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Abstract 

We discuss the potential impacts on the CMSSM of future LHC runs and possible e + e _ and 
higher-energy proton-proton colliders, considering searches for supersymmetry via Ifr events, 
precision electroweak physics, Higgs measurements and dark matter searches. We validate 
and present estimates of the physics reach for exclusion or discovery of supersymmetry via 
If T searches at the LHC, which should cover the low-mass regions of the CMSSM parameter 
space favoured in a recent global analysis. As we illustrate with a low-mass benchmark 
point, a discovery would make possible accurate LHC measurements of sparticle masses using 
the MT2 variable, which could be combined with cross-section and other measurements to 
constrain the gluino, squark and stop masses and hence the soft supersymmetry-breaking 
parameters mo, mp 2 and A 0 of the CMSSM. Slepton measurements at CLIC would enable mo 
and TO 1/2 to be determined with high precision. If supersymmetry is indeed discovered in the 
low-mass region, precision electroweak and Higgs measurements with a future circular e + e - 
collider (FCC-ee, also known as TLEP) combined with LHC measurements would provide 
tests of the CMSSM at the loop level. If supersymmetry is not discovered at the LHC, is 
likely to he somewhere along a focus-point, stop coannihilation strip or direct-channel A/H 
resonance funnel. We discuss the prospects for discovering supersymmetry along these strips 
at a future circular proton-proton collider such as FCC-hh. Illustrative benchmark points 
on these strips indicate that also in this case FCC-ee could provide tests of the CMSSM at 
the loop level. 
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1 Introduction 


The first run of the LHC at 7 and 8 TeV has framed the agenda for its future runs, and for 
possible future colliders. The CMS and ATLAS experiments have discovered a Higgs bo¬ 
son [I], but have found no sign of supersymmetry or any other physics beyond the Standard 
Model |2|,|3j. Present and future studies of the Higgs boson can be used to constrain scenar¬ 
ios for new physics, as can other high-precision low-energy measurements and cosmological 
constraints. We address in this paper the prospects for discovering supersymmetry during 
future runs of the LHC at 13/14 TeV in light of the indirect information currently provided 
by the Higgs and other measurements, and consider possible scenarios for discovering or 
measuring supersymmetry at proposed future linear and circular colliders, either directly or 
indirectly, showing how the various colliders may complement each other. 

Our study is within the minimal supersymmetric extension of the Standard Model with 
soft supersymmetry-breaking parameters constrained to be universal at a high input scale, 
the CMSSM |4,5|. This model is not imposed by top-down considerations based on string, 
M- or F-theory, nor is it required by bottom-up considerations such as limits on flavour¬ 
changing neutral interactions. However, it is the simplest supersymmetric model, so its 
phenomenology is relatively unambiguous. As such, it provides a convenient benchmark for 
considering the interplay between different high-energy colliders. 

One of the most important constraints that we take into account is the bound on the 
density of cold dark matter, which provides interesting constraints on the parameters of the 
CMSSM. In particular, requiring that the relic density of the lightest supersymmetric particle 
(LSP), assumed here to be the lightest neutralino % |6|, falls within the range allowed by 
astrophysics and cosmology can be used to provide important constraints, including upper 
limits, on the soft supersymmetry-breaking mass parameters in the CMSSM, and hence 
sparticle masses (7-91. 

The LHC measurement of the Higgs mass already provides a significant constraint on 
the parameter space of the CMSSM, favouring sparticle masses that are consistent with the 
non-observation of supersymmetric particles at the LHC in Run 1 (8-15 . The starting-point 
for our analysis is a recent global fit to the CMSSM model parameters 112), using these 
measurements as well as precision electroweak and flavour observables, as well as direct 
constraints on the interactions of the LSP with ordinary matter. 

In order to evaluate the potential of future LHC runs to probe the CMSSM, we extrap¬ 
olate the sensitivities of gluino, squark and stop searches at LHC Run 1 at 7 and 8 TeV 
to estimate LHC capabilities with 300 or 3000/fb of data at 13/14 TeV. We find that such 
data sets should permit the LHC experiments to discover supersymmetry if it has CMSSM 
parameters within the low-mass region favoured by the global fit [l2|. Assuming optimisti¬ 
cally that they are given by the best current fit in this low-mass region, we then discuss how 
accurately the LHC experiments could measure the gluino, squark and stop masses, and 
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hence the CMSSM soft supersymmetry-breaking parameters mo,mp 2 and A 0 . 

In this optimistic scenario where Nature is described by the CMSSM in the low-mass 
region, experiments at the proposed CLIC e + e _ collider at 3 TeV in the centre of mass |16,17| 
would be able to produce and measure very accurately the masses and other properties of the 
sleptons and the lighter gauginos, enabling, for example, high-precision determinations of 
the soft supersymmetry-breaking parameters mo and mp 2 of the CMSSM. An e + e _ collider 
with an energy 1 TeV could also explore parts of the low-mass region, e.g., pair-producing 
the lighter stau slepton at the low-mass best-fit point. On the other hand, e + e _ colliders 
with energies below 500 GeV in the centre of mass would not be able to produce and measure 
sparticles directly. 

As we discuss, measurements of Z-boson 18 and Higgs couplings 119-21 do not as yet 


provide strong supplementary constraints on supersymmetric models such as the CMSSM. 
However, future higher-precision measurements could be used to constrain the CMSSM pa¬ 
rameters indirectly. In particular, if Nature is indeed described by the CMSSM with pa¬ 
rameters in the low-mass region, measurements of the Z and Higgs boson at the proposed 
high-luminosity circular e + e _ collider FCC-ee (TLEP) 1221 could be used, in conjunction 
with the LHC measurements, to test this supersymmetric model at the quantum level, as 
we illustrate in the specific example of the best-fit low-mass point from 12 . As an aside, we 


also show how, again in the optimistic low-mass scenario, high-precision Z measurements at 
FCC-ee (TLEP) could be used to probe models of supersymmetric grand unification. 

On the other hand, in the pessimistic scenario where the LHC does not discover super- 
symmetry, but only establishes 95% CL lower limits on particle masses, we consider the 
prospects for discovering supersymmetry directly at a future higher-energy circular proton- 
proton collider such as FCC-hh 23,24 , or finding indirect evidence for supersymmetry via 
high-precision e + e _ measurements. Within the CMSSM, high-scale supersymmetric models 


can be found along narrow strips where stop-neutralino coannihilation is important 15 


or in the focus-point region 14 , and we analyse the prospects of direct and indirect mea¬ 
surements along these strips. Studies of illustrative benchmark points along these strips 
indicate that the combination of direct FCC-hh and indirect FCC-ee measurements could 
test supersymmetry at the loop level also in this pessimistic case. 

The layout of this paper is as follows. In Section 2 we discuss the extrapolations of 
current LHC sparticle search sensitivities to future LHC runs. Then, in Section 3 we discuss 
possible LHC measurements of particle masses in the optimistic low-mass best-fit scenario. 
Section 4 contains our discussion of e + e _ probes of supersymmetry in this optimistic scenar¬ 
ios, including direct searches at CLIC as well as indirect constraints due to high-precision Z 
and Higgs measurements at FCC-ee (TLEP). The pessimistic high-mass scenarios in which 
the LHC does not discover supersymmetry are discussed in Section 5, where we consider the 
prospects for direct discovery with FCC-hh as well as indirect measurements with FCC-ee 
(TLEP). Finally, our conclusions are summarised in Section 6. 
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2 Extrapolations of Current LHC Sparticle Search Sen¬ 
sitivities to Higher Energy and Luminosity 

The baseline for our studies is provided by a recent global fit to the parameters of the 
CMSSM [l2| D In addition to the ATLAS search for jets + I£t events with ~ 20/fb of 
8 TeV data |2,3], these global fits included the measurement of rip, 11,301 (which was related 
to the CMSSM parameters via calculations using FeynHiggs 2.10.0 |31,32|), electroweak 
precision observables and — 2 (33j, precision flavour observables including 6 —>■ sy [341 and 
B s ,d —t l l+ /r [35 -38 , and dark matter observables including the direct LUX constraint on 


dark matter scattering 39 and the total cold dark matter density |40|. These measurements 


were combined into a global y 1 2 likelihood function, whose projection on the (mo, mp 2 ) plane 
of the CMSSM is displayed in Fig. [l] In this and subsequent figures, we marginalise over 
the other CMSSM parameters tan/5 and A 0 . We display in red and blue, respectively, 
Ay 2 = 2.30 and 5.99 contours (which we use as proxies for 68% and 95% CL contours). For 
each set of (mo,ffli/ 2 ) values within these contours, there is some choice of tan/5 and A (} for 
which Ay 2 < 2.30 or 5.99, respectively, and outside these contours there are no choices of 
tan /5 and A 0 that satisfy these conditions. In the figure, a low-mass “Crimea” region and a 
high-mass “Eurasia” region can be distinguished. The former consists of points in the stau 
coannihilation region, and the latter includes points along rapid H/A annihilation funnels, 
and along the high-mass focus-point and stop coannihilation strips we discuss in Section 6. 
We also show as a filled green star a representative best-fit point in the low-mass region, 
whose parameters are fisted in Table [l] In the low-mass region, g ll: — 2 makes a significantly 
smaller contribution to the global y 2 function than in the high-mass region, although the 
CMSSM and related models could not by themselves resolve the discrepancy between the 
experimental measurement and the theoretical calculation within the Standard Model [/} We 
discuss later characteristics of points in the high-mass ‘Eurasia’ region: the y 2 likelihood 
function is relatively flat across this region, and there is no well-defined best-fit point that 
is favoured strongly with respect to other points. 

In Fig. [I] we also show as a black line the 95% CL exclusion contour in the (m 0 ,mi/ 2 ) 
plane established by ATLAS searches for jets +^t events with ~ 20/fb of data at 8 TeV |2|. 
This exclusion was derived within the CMSSM with tan /5 = 30 and (in our sign convention) 
A 0 = 2 mo, but studies have shown that the limit is relatively insensitive to the values of 
tan/5 and A 0 11 . The ATLAS 95% CLs contour intersects the 95% CL contours found in 


the global fit, reflecting the importance of other observables in the global fit. For example, as 


1 This paper also contains a global fit to the NUHM1 25,26 , and a fit to the NUHM2 [26j[27| can be 

A further update of the 


found in 128 


together with minor updates of these CMSSM and NUHM1 fits. 

These updates do not impact qualitatively the analyses presented in 


CMSSM analysis can be found in 29 
this paper. 

2 This discrepancy can be resolved in a model that relaxes the assumption of GUT-scale universality in 


the soft supersymmetry-breaking parameters 29 . 
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Figure 1: The (mo, mi/ 2 ) plane in the CMSSM. The Ay; 2 = 2.30 (68% CL) and 5.99 (95% 
CL) regions found in recent global fits are bounded by solid red and blue lines, respectively. 
The best-fit point in the low-mass ‘Crimea’ regions is indicated by a filled green star. Also 
shown as solid black (purple, green) lines are the sensitivities of LHC Ifix searches for ex¬ 
clusions at the 95% CLs with 20/fb of data at 8 TeV (300, 3000/fb of data at If TeV). The 
purple contour is expected to coincide (within uncertainties) with the 5-a discovery contour 
at the LHC with 3000/fb of data at If TeV. 


Model 

Location 

m 0 

(GeV) 

mi/2 

(GeV) 

Ao 

(GeV) 

tan /3 

CMSSM 

Low-mass 

670 

1040 

3440 

21 


Table 1: Representative low-mass best-fit point found in a recent global CMSSM fit [121, using 
the ATLAS jets + Ifr constraint & and the combination of the LHCb and CMS p*7^ 
constraints on B S: d —>■ p + p~ [381, and using FeynHiggs 2.10.0 [321 to calculate m / l . 


already mentioned, — 2 tends to favour relatively low values of mo and in\/ 2 - On the other 
hand, the measurement of tends to favour values of m 0 and rrii / 2 beyond the ATLAS 
contour. 

We use a simple procedure to estimate the sensitivities of future collider searches ex- 
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ploiting th efa signature accompanied by jets (possibly 6-tagged) and/or leptons at higher 
centre-of-mass energies and luminosities. We scale the 95% CL exclusion or 5-cr discovery 
contours of the searches at 8 TeV to different luminosity and energy scenarios by assuming 
that the signal efficiency and background suppression of the current 8- Te V searches remain 
unchanged. Maintaining the present performance of the searches is motivated by the AT¬ 
LAS and CMS upgrade programmes, and is defined by both experiments as one of the main 
upgrade goals. The assumption was also used in several studies for Snowmass and ECFA 
(see e.g. |41,42 ) as well as to project collider limits for Dark Matter searches (43,44 . It also 


forms the basis of the Collider Reach (45] tool, which reports dedicated studies showing good 
agreement between this extrapolation approach and results obtained from a full simulation. 

We caution, however, that various effects could invalidate our assumption. For exam¬ 
ple, the signal-to-background ratio could vary with the centre-of-mass energy and with the 
number of pile-up events, which is correlated with the luminosity. Indeed, our extrapolation 
of the LHC sensitivities with 300/fb and 3000/fb of integrated luminosity at 14 TeV in the 
centre of mass is somewhat less conservative than ATLAS estimates of their exclusion sensi- 
However, we have been able to verify that our simple assumption gives similar 


tivities 


46 


results to Snowmass estimates of the possible sensitivities of higher-energy colliders based on 
simplified model searches at the LHC with ~ 20/fb of data at 8 TeV (44], and we consider 
our assumption a reasonable objective for future experimental analyses to target. 

We consider in this paper the following LHC sparticle searches: searches for events with 
jets and missing transverse energy, fa, possibly accompanied by leptons and with some 
jets 6-tagged, dedicated searches for light stop squarks t —> x + c , an d monojet searches. 
Using our simple assumption for a number of current LHC searches, we calculate cross- 
sections at higher LHC centre-of-mass energies with Pythia 8 [ 47,48], using as default the 
MSTW2008NLO parton distribution functions |49 We then require that the products of 
the integrated luminosity with the cross-section be the same as for the 8 TeV LHC data. In 
this way, we extrapolate current LHC 95% CLs exclusion limits to higher LHC energies and 
luminosities, as well as possible future colliders with 3000/fb at 33 and 100 TeV, as seen in 

Table 00 

Fig. m displays as purple and green lines, respectively, our extrapolations within the 
CMSSM of the current ATLAS 95% CLs limit from searches for jets + fa events with 
~ 20/fb of data at 8 TeV to LHC searches at 14 TeV (LHC 14) with 300/fb and 3000/fb of 
integrated luminosity. (We note that the ATLAS study (46] found that the 5-cr discovery 
contour for 3000/fb almost coincides with the 95% CLs exclusion contour for 300/fb.) Within 
the CMSSM, the ATLAS search for jets + fa events is the most sensitive for m 0 /mi/ 2 < 


3 We have verified that very similar results can be obtained using the NNPDF2.3LO, option 14 parton 
distribution functions [50 . 

4 There are also extrapolations available for slepton and chargino searches. However, in those cases the 
interpretations of the searches are more delicate 


29 


and we do not discuss them here. 
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LHC 

HE-LHC 

FCC-hh 

Search 

signature 

8 TeV 
20/fb 

14 TeV 
300/fb 

14 TeV 
3000/fb 

33 TeV 
3000/fb 

100 TeV 
3000/fb 

G 9 bb X) 2 ( m g) 

1300 

2540 

2990 

6080 

14700 

it* (m t ~) 

650 

1350 

1740 

3260 

7020 

(t -> ex) 2 (mf) 

240 

530 

780 

1320 

2510 

CMSSM 

(m 0 ,m 1/2 ) 

ra 0 = m 1/2 

(800, 800) 

(1610, 1610) 

(1860,1860) 

(4080, 4080) 

(10800, 10800) 

mo = 2.5m 1 / 2 

(1500, 600) 

(2950, 1180) 

(3390, 1360) 

(7310, 2930) 

(19000, 7600) 


Table 2: Extrapolations of current LHC searches with ~ 20 /fb of luminosity at 8 TeV 
to higher energies and luminosities, assuming sensitivities to the same numbers of signal 
events. The first five rows of the Table are possible 95% CL exclusion sensitivities derived 
from searches for specific sparticle pair-production processes, as indicated, and the numbers 
correspond to the sparticle masses in GeV. The last two rows are for rays in the (m 0) m i/ 2 ) 
plane, as indicated, and the numbers correspond to the possible 95% CL exclusion limits on 
mo and mi/ 2 - 


2, with other searches becoming more important at larger mo/mi/ 2 - We return later to 
extrapolations of mono jet searches and dedicated searches for light stop squarks, which are 
important for our studies of FCC-hh. 

We see that the low-mass ‘Crimea’ region lies within the purple (95% CLs exclusion with 
300/fb or 5-cr discovery with 3000/fb at 14 TeV) contour where mo < whereas the high- 
mass ‘Eurasia’ region lies largely beyond the purple contour. Based on these comparisons 
between the extrapolated LHC sensitivity and current fits within the CMSSM, we have 
chosen for further study two scenarios for the outcome of the LHC searches with 3000/fb at 
14 TeV. 

• An ‘optimistic’ scenario in which the LHC discovers supersymmetry in the ‘Crimea’ 
region, and for definiteness we assume that its parameters coincide with those at the 
the representative low-mass best-fit point in Table [lj 

• A ‘pessimistic’ scenario in which the LHC discovers no evidence for supersymmetry, in 
which case the supersymmetry-breaking parameters must he somewhere in ‘Eurasia’. 

The following sections contain discussions of the interplay between the various colliders 
in these scenarios. 
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3 LHC Measurements of Supersymmetry in the 
Optimistic Scenario 

Assuming that Nature is described by supersymmetry at the CMSSM low-mass best-fit 
point, the sparticle mass spectrum is determined, as illustrated in Fig. [2| The most relevant 
sparticles for searches at the LHC are those with the highest production cross sections, 
namely squarks and gluinos. At the best-fit point, the mass of a generic right-handed u, d, s, c 
or b squark is calculated to be rriq R ~ 2080 GeV, and the lighter stop squark has a mass 
mq ~ 1020 GeV. Also, nig ~ 2280 GeV and the lightest neutralino mass rn x ~ 450 GeV. 
The lighter stau mass m fl is only very slightly heavier: at the best-fit point and the rest of 
the low-mass region stau-y coannihilation is responsible for bringing the relic density into the 
range allowed by cosmology. In the following we consider the possible LHC measurements of 
generic right-handed squarks qn, gluinos g and the lighter stop A, with either 300 or 3000/fb 
of luminosity at LHC 14. We assume that experiments at the LHC discover supersymmetry 
with the mass spectrum characteristic of the best-fit point shown in Fig. [2j and ask how 
accurately its parameters can be measured. 



Figure 2: The spectrum at the best-fit point in the CMSSM 121, whose parameters are listed 
in Table [!} The magnitudes of the branching ratios for sparticle decays into different final- 
state particles are represented by the strengths of the dashed lines connecting them. 
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3.1 Gluinos and Squarks 


We estimate first the potential resolution with which the gluino and squark masses could 
be measured. For this purpose, we consider three contributions to the determination of 
these model parameters: measurements of the total cross-section, the distribution in the 
MT2 variable 51 , and the spectator jet energies in g —y q + q R decay. Fig. [3] shows how 
the total cross-section for strongly-interacting particle production at LHC 14 obtained from 
Pythia depends on the gluino mass trig (left panel) and the squark mass rriq R (right panel), 
expressed as functions of the mass differences AM relative to the low-mass best-fit values 
in the CMSSM [^J We see that the dependence of the cross-section on nig is much weaker 
than that on rriq R . In the following we combine the information that can be derived the 
cross-section with that obtainable from an analysis using the MT2 variable. 




AM/GeV AM/GeV 

Figure 3: The sensitivities of the total sparticle cross-section to nig (left panel) and nriq R 
(right panel), expressed as functions of the mass differences AM relative to the low-mass 
best-fit values in the CMSSM. 

In order to assess how MT2 measurements could contribute to constraining the gluino 
and squark masses, we construct a set of MT2 templates for various values of these and 
the neutralino mass, and fit these templates to a simulation of the prospective MT2 dis¬ 
tribution for the central best-fit values of the masses. For this analysis, we first matched 
the reconstructed jets from the Pythia output to the squarks, gluinos and neutralinos at 
the generator level. We then applied the same logic to construct MT2 as in experimental 
papers, treating the neutralinos as Ifr and the decay products of the squarks and gluinos 
as the jets. Thus, this approach does not consider combinatoric effects as could arise in a 
study that used a full detector simulation. We use the same kinematic specifications for the 

5 We note that the relative fractions of the gluon-gluon, gluon-squark and squark-squark final states vary 
continously in these plots. 







search regions as in the published 8 TeV search |2|, and assume that the sensitivity remains 
the same for 14 TeV 0 

Fig. [4] displays prospective histograms of the MT2 distributions obtained from simulations 
using Pythia 8 |47[|48] and the MSTW2008NLO parton distribution functions j49| for 
different values of rrig (upper panel), the right-handed squark mass nig (lower panel) Q In 
both cases, we compare the distribution for the nominal mass at the best-fit point with 
the corresponding distributions for values of the mass deviating from the nominal value by 
±300 GeV, keeping the other sparticle masses fixed. In the gluino case, we see that the 
MT2 histogram for the nominal value nig = 2280 GeV (in red) is very similar to that for 
the —300 GeV choice (in blue), whereas the histogram for the ±300 GeV choice is less 
similar. The reverse is true for the squark case (middle panel): here the nominal histogram 
for rriq R ~ 2080 GeV (red) is more similar to that for the ±300 GeV choice (green), and less 
similar to that for the —300 GeV case (blue). 

The plots in Fig. [4] were obtained by recalculating the full Pythia output as the squark 
and gluino masses were varied around the best-fit CMSSM point. In some cases, the variation 
changed the ordering of the squark and gluino masses, leading to substantial changes in the 
MT2 distribution, e.g., in the AM = —300 GeV case in the upper panel of Fig. [4] (green 
histogram) and in the AM = ±300 GeV case in the lower panel (blue histogram). The 
changes in the shapes of the MT2 distributions were less important when the mass ordering 
stayed the same. In addition, the variations in the shape of the MT2 distribution include 
the effects of changes in the relative production rates of gg, gq , qq and qq final states arising 
from the mass variation. 

In order to estimate the uncertainties in measurements of sparticle masses that could 
be possible at the LHC, we have performed fits to the simulated data for varying amounts 
of integrated luminosity. We use these to estimate the 68% CL ranges of mass estimates 
obtainable with either 300 or 3000/fb of integrated luminosity. As seen in Fig.[3J the changes 
in the MT2 distributions for gluino and squark mass changes of ±300 GeV are quite different, 
so we do not expect symmetric Gaussian uncertainties, and we note that the same is true for 
the projected cross-section measurements shown in Fig. [3} Combining these with the MT2 
measurements, we find the ;y 2 distributions as functions of nig and rrig R shown in Fig. [5] in 
the left and right panels, respectively. The y 2 functions are evaluated as 


X' 


\m) = V 

n ^ 


2=1 


(Mm) - 
af 


(i) 


where rh is the nominal mass, the Ni are numbers of events in the simulation and cq is the 

6 We consider this to be a conservative assumption, as the signal-to-background ratio is likely to improve. 

7 Here and subsequently, we include in our simulations the Standard Model backgrounds from Z°, 
it and single-£ production. However, we do not embark on a simulation of either ATLAS or CMS, since the 
purpose of our exploratory study is to give a first feeling for what might be possible in future LHC runs, 
and suitable detector simulations for the high-luminosity LHC are not available. 
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Figure 4: Simulations for lf-TeV collisions, using Pythia 8 [47,48] and including Standard 
Model backgrounds, of the distributions in the MT2 variable for (upper panel) the nominal 
value of the gluino mass at the low-mass CMSSM best-fit point, nig ~ 2280 GeV (blue 
histogram), and gluino masses differing by ±300 GeV (green and blue histograms), and 
similarly for (lower panel) the nominal value of the squark mass nig ~ 2080 GeV and values 
±300 Ge V. In both cases, we fix the other sparticle masses to their nominal best-fit values, 
assuming in particular that the LSP mass m x = 450 GeV. The inserts show the integrated 
luminosities at 14 TeV that would be required to distinguish at the 3-a level between the best 
fit and other models with the indicated mass shifts AM. 

statistical error in each bin for the assumed luminosity, and the sum over i = 1, ...n includes 
all the bins in the histograms added in quadrature. The upper row of panels is for 300/fb of 
integrated luminosity, and the lower row is for 3000/fb of integrated luminosity. On the basis 
of this analysis, we estimate the following fit uncertainties with 300/fb of data at 14 TeV: 

300/fb: A nig = (-270, ± ...) GeV , 

A m~ qR = (-100, ±110) GeV. (2) 
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where the ... indicate that these measurements provide no useful upper limit on rrig , and 
with 3000/fb: 


3000/fb : Arrig = (-110,+150) GeV, 

A rrig R = (-30, +35) GeV. (3) 

These uncertainties do not include a potential systematic effect from jet energy scale uncer¬ 
tainties. However, as we expect these to be at the level of 10% or below, their overall impact 
is expected to be sub-dominant. 




A M/GeV 


AM/GeV 




A M/GeV 


AM/GeV 


Figure 5: The % 2 functions for nig (left panels) and niq R (right panels), as estimated from 
cross-section and MT2 measurements with 300/fb (upper panels) and 3000/fb (lower panels). 


The upper and lower panels of Fig. [4] show that the mass difference nig — m,q R is poorly 
constrained by the MT2 measurement, and this is reflected in the asymmetric \ 2 distribu¬ 
tions seen in Fig. [5| However, there are many other possible measurements at the LHC. 
In particular, we have considered the extra information that could be obtained from mea¬ 
surements of the (relatively) soft jet emitted in the decay g —» qR + q, which would be 
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p T in g -» q R + q decay 


monochromatic in the gluino rest frame. Fig. [6] displays scatter plots of the Pt of the jet 
emitted in —> q + % decay (horizontal axis) and the jet emitted in g —» + q decay 
(vertical axis) based on simulations of g pair-production. The left panel is for the best-fit 
values of nig and m rjR , and the right panel is for the same value of nig but with nriq R reduced 
by 300 GeV: the plots are clearly distinct. 



Figure 6: Scatter plots of the pr (in GeV) of the jet emitted in qR —» q + x decay (horizontal 
axis) and the jet emitted in g —* qR + q decay (vertical axis) resulting from a simulation of 
gluino pair-production at the LHC at If TeV. Left panel: for the best-jit g and qn masses. 
Right panel: for the best-jit values of mg but with rn,q R reduced by 300 Ge V. 


Fig. [7] displays the spectrum of the ‘soft’ jet in the same two cases: in the left panel with 
the best-fit g and qp, and in the right panel with the same value of trig but with eri ( j R reduced 
by 300 GeV. These can clearly be distinguished with a high degree of confidence. We do not 
display the corresponding distribution with m g - fi increased by 300 GeV, since in this case the 
qR is heavier than the g and there is no ‘monochromatic’ supplementary jet in gluino decay. 

We assume that the jet energy in g —>■ qR + q decay can be measured with an accuracy of 
50 GeV. This information can then be combined with the cross-section and MT2 distribution 
discussed earlier to estimate 68 and 95% CL regions in the ( m,q R ,mg ) plane. These are shown 
shaded pink and blue, respectively, in Fig. [8] for 300/fb of integrated luminosity (upper left 
panel) and for 3000/fb of integrated luminosity (upper right panel). As in Fig. [lj the 
low-mass portions of the solid red and blue contours outline the Crimea region and the high- 
mass portions correspond to the Eurasia region Finally, the solid [dashed] magenta lines 
(darker and lighter) show the 5-cr discovery (95% CL exclusion) reaches of the LHC with 
300 (3000)/fb. The lower panels of Fig. [8] show as solid red and blue lines the 68 and 95% 
CL contours from fits combining the prospective LHC measurements with the recent global 
fit 1121 (whose CL contours are displayed as dashed lines in these panels). 

These prospective measurements can be projected onto the (mo, mi/ 2 ) plane of the 
CMSSM, as seen in the upper panels of Fig. [9j also for 300/fb of integrated luminosity 
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Figure 7: Simulations of of the distributions of the quark pr (in GeV) from g pair-production 
at the LHC at If Te V followed by g -4 <)r + q decays. Left panel: for the best-fit g and qR 
masses, and (right panel) for the same value of nig but with mq R reduced by 300 GeV. 


(upper left panel) and 3000/fb of integrated luminosity (upper right panel). Formally, the 
corresponding numerical 68% CL uncertainties are 


and 


300/fb : A m 0 = 
3000/fb : A m 0 = 


300/fb : Amp2 = 
3000/fb : Amp2 = 


(-670, +620) GeV, 
(-670, +220) GeV, 


(-140,+100) GeV, 
(-90, +20) GeV. 


( 4 ) 


( 5 ) 


The upper panels of Fig. [9] also show that these numbers imply non-trivial correlations be¬ 
tween m 0 and m ,/ 2 • They also show, as solid red and blue lines, respectively, the boundaries 
of the 68 and 95% CL regions found in the recent global analysis of current data [ 121. We 
see that the prospective future LHC measurements could provide information that would be 
complementary to this global fit 110 

The lower panels of Fig. [9] show the results (solid red and blue lines) of combining these 
LHC measurements with the recent global fit [12] (dashed red and blue lines). We see that 
the LHC measurements would reduce substantially the sizes of the 68 and 95% CL regions 
already with 300/fb, and that the prospective 3000/fb measurements would be particularly 
powerful in this regard. 

8 The main effect of the measurement of the supplementary jet in g —>• qyi + Q decay is to truncate the 
preferred strip in the 300/fb case: it has no visible effect in the 3000/fb case. 


13 














Figure 8: The upper panels show the 68 and 95% CL regions (shaded pink and blue, re¬ 
spectively) in the (rriq R ,mg) planes obtained from cross-section, MT2 and supplementary jet 
measurements at LHClf with 300/fb (left panel) and 3000/fb (right panel). These regions 
are superposed on the best-fit point (green star) and the 68 and 95% CL regions found in a 
recent global fit to the CMSSM (12) (solid red and blue lines), and the magenta lines show 
the prospective capabilities of the LHC to exclude at the 95% CL (dashed) supersymmetry 
or discover it at the 5-a level (solid) with 300/fb or 3000/fb (darker and lighter lines). In 
the lower panels we show as solid red and blue lines the results of fits combining the LHC 
measurements with the recent global fit (here shown as dashed lines). 


3.2 Stop Measurements 


We have also considered the possible accuracy in measuring via t\ + t\ production at 
the LHC at 14 TeV. The left panel of Fig. [TO] shows the sensitivity of the total stop pair- 
production cross-section to : we see that over the displayed range it is greater than those 
to the ntg and niq R , that were shown in Fig. [3j The right panel of Fig. 10 shows histograms 


of MT2 for the nominal mass ~ 


1020 GeV at the representative low-mass best-fit point 
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Figure 9: The upper panels show the 68 and 95% CL regions (shaded pink and blue, re¬ 
spectively) in the (mo, 7711 / 2 ) planes obtained from, cross-section, MT2 and supplementary jet 
measurements at LHClf with 300/fb (left panel) and 3000/fb (right panel). These regions 
are superposed on the 68 and 95% CL regions found in a recent global fit to the CMSSM [12] 
(red and blue lines). In the lower panels we show as solid lines the results of fits combining 
the LHC measurements with this global fit (here shown as dashed lines). 


and for choices differing by ±300 GeV. These cases are quite distinct, as is also seen in the 
inset, which displays the luminosities required for 3-cr discrimination between the nominal 
value of mq and selected larger or smaller values. 

The left panel of Fig. [Tl] displays the shape of the unit-normalised tt invariant mass 


distribution resulting from a simulation of such events using Pythia 8 47, 48 and the 


MSTW2008NLO parton distribution functions 49 , produced with the nominal CMSSM 


best-fit values of m/ i = 1020 GeV and nig = 2280 GeV (green histogram), compared with 
the Standard Model background (black histogram), which is sharply peaked at low invariant 
masses close to the tt threshold. Also shown in Fig. [lT] are the invariant-mass distributions 
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Figure 10: Left panel: The sensitivity of the total stop pair-production cross-section to m k . 
Right panel: Simulations for lf-TeV collisions of the distributions in the MT2 variable for 
the nominal value of the lighter stop mass m k = 1020 GeV and values ±300 GeV, with the 
other sparticle masses fixed to their nominal best-fit values. The insert shows the integrated 
luminosities at 14 TeV that would be required to distinguish at the 3-a level between the best 
fit and other models with the indicated mass shifts AM relative to the value at the low-mass 
best-fit point. 


for g masses 300 GeV above (red histogram) and 300 GeV below (blue histogram) the 
nominal value of m k . As expected, the higher (lower) mass gives a longer (shorter) tail in 
the invariant-mass distribution. On the other hand, as we see in the right panel of Fig. 0 
that the invariant i\ ± t\ mass distribution in g decays is almost independent of rn k for fixed 

m-g. 


Combining the cross-section, MT2 and tt invariant-mass measurements, we find the ;y 2 
distributions as functions of m k shown in Fig. 12. The left panel is for 300/fb of integrated 
luminosity, and the right panel is for 3000/fb of integrated luminosity. We find the following 
fit uncertainties with 300/fb or 3000/fb of data at 14 TeV: 


300/fb: A m h = (-30, ±50) GeV, 

3000/fb: Am k = (-10, ±15) GeV. (6) 

As in the previous cases, these uncertainties should be convoluted with a systematic jet 
energy scale uncertainty of ~ 10%. 

The uncertainties (|6j) may be used to estimate the corresponding uncertainties in the 
trilinear soft supersymmetry-breaking parameter A 0 , by comparing the stop mass (which is 
very sensitive to A (l ) with the squark and gluino masses (which are insensitive to A 0 ). The 
effect of marginalising over the latter masses can be incorporated by assuming that rn 0 and 
mi /2 have their best-fit values, as is also the case for tan/5. In estimating the uncertainty 

in A 0 , we incorporate the correlation between A 0 and // that is imposed by the electroweak 
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Distribution (arbitrary units) 



Figure 11: The unit-normalised, tt invariant-mass distribution resulting from a simulation 
of t\ + ti production at the LHC at If TeV. Left panel: for the best-fit g and t\ masses of 
2280 and 1020 GeV (green histogram), compared with the Standard Model background (black 
histogram) and simulations with g masses 300 GeV above (red histogram) and 300 GeV below 
(blue histogram) the nominal value of nig. Right panel: similarly for the best-fit g and t\ 
masses (green histogram), compared with the Standard Model background (black histogram) 
and simulations with t\ masses 300 GeV above (red histogram) and 300 GeV below (blue 
histogram) the nominal value of . 



Figure 12: The x 2 functions for rn^ estimated from LHClf measurements with 300/fb (left 
panel) and 3000/fb (right panel). 


vacuum conditions within the CMSSM, finding 

300/fb: AA 0 = (+80,-150) GeV, 

3000/fb: AA 0 = (+30,-40) GeV. (7) 

We emphasise again that these uncertainties do not take into account the jet energy scale 
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uncertainty, which we expect to be subdominant. 

As a final point in this Section, we comment on the magnitudes of some of the branching 
ratios for sparticle decays that are represented by dashed lines in Fig. [2} The earlier analysis 
of g —* q + q decays exploits the fact that this decay mode is dominant if rrig > m,q R , 
as reflected in the boldness of the dashed line connecting the g and <%> states. We draw 
attention to the decay fa ti + h , which is also dominant, having a branching ratio of 77% 
represented also by a bold dashed line. This implies that about 50% of t 2 t 2 events would 
contain, in addition to a tt pair, a pair of high-p r Higgs bosons and substantial missing 
transverse energy. Typical boost factors for the Higgs bosons would be ~ 5. A detailed 
exploration of this experimental signature lies beyond the scope of this paper. 


4 e + e Probes of Supersymmetry in the Optimistic Sce¬ 
nario 

In the low-mass ‘optimistic’ CMSSM scenario there would be interesting opportunities for 
both direct and indirect precision probes of supersymmetry at an e + e _ collider, which we 
now explore. 


4.1 Direct Sparticle-Pair Production 


The most direct possibility would be pair-production and measurement of electroweakly- 
interacting sparticles. The next-to-lightest supersymmetric particle (NLSP) is expected, in 
generic regions of the CMSSM parameter space, to be the lighter stau slepton f\. Accordingly, 


Fig. 13 displays, superimposed on the same CMSSM (m 0 , mi/ 2 ) plane discussed previously, 
contours showing where it is possible at the 95% CL to attain = 500 GeV (green), 
the largest mass that could be pair-produced with an Ecm = 1 TeV linear collider, and 
1500 GeV (black), the largest mass that could be pair-produced with an Ecm = 3-TeV 
linear collider such as CLIC. These contours are restricted to the regions within the 68 and 
95% CL regions found in the recent global fit |12|, where the CMSSM parameter space is 
well sampled. We see that the = 500 GeV line crosses the ‘Crimea’ region, whereas 

the m fl = 1500 GeV lines reach deep into the ‘Eurasia’ region. In particular, the low-mass 
best-fit point in the CMSSM lies within the = 500 GeV reach of a 1-TeV e + e _ collider. 
In the low-mass ‘Crimea’ region, the cold dark matter density is brought into the range 
acceptable to cosmology by coannihilation with the stau, so the m fl < 500 GeV contour 
has rrii /2 almost constant. On the other hand, in the ‘Eurasia’ region other mechanisms 
such as neutralino annihilation via direct-channel heavy Higgs poles come into play, and the 
m fi < 1500 GeV contour has a more complicated shape. 

Within the specific CMSSM model studied, a 500-GeV e + e _ collider would very likely 
be unable to detect directly any supersymmetric particles. This is because the contour for 
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Figure 13: Contours where it is possible to attain at the 95% CL rrif l = 500(1500) GeV, 
indicated by solid green (black) lines, are overlaid on the (mo^i/ 2 ) plane in the CMSSM, 
with the same CL contours and best-fit point from a global fit as displayed previously in 
Fig- 0 


rftf, = 250 GeV, the largest mass that could be pair-produced with an Eqm = 500-GeV 
linear collider, would lie at mi/ 2 < 600 GeV, which is outside the 95% CL contour in the 
(mo, m 1 / 2 ) plane shown in Fig. 13 A similar conclusion could be drawn from the lower right 


panels of Figs. 5 and 13 of 1121, where we see that Ay 2 > 9 for rn^ < 250 GeV^J 

If slepton-antislepton pair-production is accessible at an e + e _ collider, many very precise 
direct measurements become possible. Two benchmark supersymmetric scenarios were anal¬ 
ysed in 1161, and the prospective accuracies for sparticle mass measurements were assessed. 
In one of these scenarios (PI), the slepton mass spectrum was very similar to that at the 
low-mass CMSSM best-fit point (see Table 1 and Fig. [2j), with masses between 1000 and 
1100 GeV. The xf and Xa masses in scenario PI were somewhat lower than in the low-mass 
best-fit CMSSM spectrum, whereas the x® mass was again very similar. Based on simula¬ 
tions of 2/ab of CLIC data at 3 TeV, the following uncertainties in sparticle masses were 
estimated: 


A mi R = 2.9GeV, Am x o = 4.6GeV, A m-± = 3.6GeV. 


"xi — — • • (8) 

The measurement uncertainty in m x 0 can be converted directly into the corresponding un¬ 
certainty in mi/ 2 : 

Ami /2 = 11 GeV. (9) 


On the other hand, this would be possible within the pMSSMIO analysis discussed in 29 . 
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Combining this uncertainty with the uncertainty in the rn f : lt measurement, one finds 

Amn = 4 GeV. 


( 10 ) 


As discussed in 16 and 17 , many precision measurements of supersymmetric particle 


masses and other properties would be possible at CLIC point, including tests of the uni¬ 
versality hypotheses of the CMSSM. However, it is already clear that CLIC could provide 
exceptional precision in the determination of CMSSM model parameters, if Nature is de¬ 
scribed by a model in the Crimea region. Moreover, the comparison between the CLIC 
determinations of the CMSSM parameters with those from the LHC discussed earlier would 
enable non-trivial checks to be made of the consistency of the CMSSM assumption of uni¬ 
versal input soft supersymmetry-breaking parameters. 


4.2 Electroweak Precision Observables 

It is also possible to obtain indirect information about supersymmetric models from elec¬ 
troweak precision observables (EWPOs), similar in principle to the information about rn t and 
rriH that were obtained previously from precision measurements at LEP and the SLC |18|. 
The left panel of Fig. [T4] displays as blue points with error bars the central values and 1-cr 


uncertainties of several such observables, as calculated in a recent global fit 52 , compared 


with their values and current individual experimental uncertainties in the Standard Model. 
Also shown (without theoretical uncertainties) are the values of these observables calculated 
at the representative low- and high-mass best-fit points in the CMSSM found in 12 . As 


is apparent from the left panel of Fig. [14] and the upper left panel of Fig. [15] the current 
experimental error in the measurement of T z is too large to provide much information about 
supersymmetric model parameters. The entire region of the CMSSM (mo,mi/ 2 ) plane cur¬ 
rently allowed at the 95% CL according to the global fit 112 is compatible with the current 
measurement of T z at the 1-cr level 18 . However, also shown in the left panel of Fig. fl4] as 


turquoise bars, are the prospective experimental errors in measurements at FCC-ee (TLEP) 
(neglecting theoretical uncertainties) 1221, normalized relative to the current experimental 
errors. It is clear that, for T z and many other electroweak precision observables, the prospec¬ 
tive FCC-ee (TLEP) uncertainties are sufficiently small to be very sensitive to deviations 
from their Standard Model values and capable of constraining supersymmetric scenarios. 

The right panel of Fig. [14] makes a similar point for precision Higgs observables, by 
comparing the estimated precisions of measurements at the LHC, the ILC and FCC-ee 
(TLEP) |22,53| (shown as colour-coded horizontal bars) with the deviations of the observ¬ 
ables from their Standard Model values that are calculated for the low- and high-mass 
CMSSM best-fit points |12|. It is clear that FCC-ee (TLEP) has the best ability to distin¬ 
guish these models from the Standard Model, as we discuss in more detail later. 

As a first example of the possible utility of the precision electroweak measurements pos¬ 
sible with FCC-ee (TLEP), we consider the optimistic scenario in which supersymmetry is 
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Figure 14: The left panel compares current measurements of electroweak precision observables 
(EWPOs) taken from a Standard Model fit / 52 / (blue, with error bars), predictions at low- and 
high-mass best-fit points in the CMSSM m (red and purple symbols) and prospective FCC- 
ee (TLEP) experimental errors fH^ / (turquoise bars). The right panel compares prospective 
measurements of Higgs branching ratios at future colliders, low- and high-mass CMSSM 
predictions (red and purple symbols) and the current uncertainties within the Standard Model 
(turquoise bars). 


within the LHC discovery range and assume, for definiteness, that the model parameters 
correspond to the current best-fit low-mass point in the CMSSM. We see in the upper left 
panel of Fig. [15] that, as already commented, the current experimental uncertainty in T z, 


namely AYz = 2.3 MeV 18 , is too large to provide significant information about CMSSM 
model parameters within the 95% CL regions displayed in Fig. [lj For this reason, all of the 
95% CL region in the upper left panel of Fig. [T|is shaded green, since it lies within one current 
standard deviation of the present measurement. On the other hand, we see in the upper right 
panel of Fig. 15 that the prospective experimental uncertainty at FCC-ee (TLEP), namely 
AY Z = 0.1 MeV 1 22 , is far smaller than the variation in Tz across even the CMSSM 68 % 
CL Crimea region. For this reason, much of the 68 % and 95% CL regions in this panel of 
Fig.0 are unshaded, since they lie more than three current standard deviations away from 
the prospective measurement. The same holds for other electroweak precision observables 
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Figure 15: The present measurement of Tz (upper left panel) (18^, and prospective FCC- 
ee (TLEP) measurements fH| / of T z (upper right), M w (lower left) and Re (lower right) 
are superposed on the preferred region of the (m 0 ,mi/ 2 ) plane in the CMSSM [12] shown 
previously in Fig. [7| The colours represent deviations from the present central value in units 
of the present LHC experimental error (upper left panel), and the deviations from the values 
at the low-mass best-fit point of the values at other points in the (m 0 , mi/ 2 ) plane in units of 
the estimated future FCC-ee (TLEP) experimental errors (other panels). 


such as Mw (prospective experimental uncertainty 0.5 MeV |22j, lower left panel of Fig. 15), 
Re (prospective experimental uncertainty 5 x 10 -5 |22|, lower right panel of Fig. 15) and 
others not shown. 

We have made a crude estimate of the impact on the recent global fit to the CMSSM 
parameters of these FCC-ee (TLEP) electroweak measurements, neglecting the inevitable 
improvements in flavour, dark matter and Higgs observables, and setting aside the direct 
measurements of sparticle masses possible at the LHC following discovery in this optimistic 
scenario. As we see in Fig. [l 6 j the electroweak precision measurements would, by themselves, 
provide very tight constraints on the CMSSM parameters m 0 and m^. 

After inclusion of the FCC-ee (TLEP) measurements 22 in this optimistic scenario, 
only a small part of the low-mass ‘Crimea’ region is allowed at the 68 or 95% CL, as seen 
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Figure 16: The prospective Ax 2 = 2.30 (68% CL) and Ax 2 = 5.99 (95% CL) contours (solid 
red and blue lines, respectively) in the (mo, m^) plane for the CMSSM (with the present 68 
and 95% CL contours shown as dashed red and blue lines, respectively), assuming that the 
electroweak precision observables are measured at FCC-ee (TLEP) to have the same central 
values as at the current low-mass CMSSM best-fit point [12] (shown as the filled green star), 
and neglecting inevitable improvements in other constraints on the supersymmetric models. 
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in Fig. 16 The impact of the FCC-ee (TLEP) measurements may be translated into the 
one-dimensional likelihood functions for various sparticle masses, shown as solid red lines in 
Fig. 17 We see that nig , rriq, rrif ] and rn^ could be estimated with interesting accuracy 
on the basis of FCC-ee (TLEP): 


nig G (1680,2480) GeV, 
nig G (1680,2280) GeV, 
m fl G (340,500) GeV, 

m~ tl G (810,1110) GeV, (11) 


whereas the nominal values at the best-fit point are 2280, 2080, 450 and 1020 GeV, respec¬ 
tively. Since, in this optimistic scenario, squarks and gluinos would have been discovered 
previously at the LHC, measurements of their masses could be compared with the estimates 
based on the FCC-ee (TLEP) measurements. Agreement would constitute a non-trivial 
test of the CMSSM at the loop level, analogous to the tests of the Standard Model made 
possible by measurements of m t and m# and their consistency with predictions based on 
LEP and SLC data |18|. Conversely, any disagreement could be interpreted as a possible 
deviation from the CMSSM assumptions of universality for the soft supersymmetry-breaking 
parameters. 


4.3 Precision Higgs Observables 


We have made a similar estimate of the potential impact of the high-precision Higgs mea¬ 


surements possible with FCC-ee (TLEP) 122], as illustrated in the right panel of Fig. 14 In 
the upper left panel of Fig. [18] we display the deviation of the present experimental value 
of BR (H —>• ZZ) from the values calculated at points within the 68 and 95% CL regions 
in the (m 0 , mi/ 2 ) plane of the CMSSM, in units of the present experimental error. In the 
other panels of Fig. [l8]we show the numbers of FCC-ee (TLEP) cr’s by which the values of 
BR (H —> ZZ) (upper right), BR (H -G WW) (lower left) and BR (H —> gg) (lower right) 
calculated at other points in the CMSSM (mo,mi/ 2 ) plane differ from the values at the 
low-mass CMSSM best-fit point. 

As in the case of the electroweak precision measurements shown in Fig. [L5j we see in the 
upper left panel of Fig. 18 that the entire 68 and 95% CL regions of the CMSSM (mo, mi/ 2 ) 


plane he within a single current LHC a of the present central value of BR (H -G ZZ), and 
hence are shaded green, and we have checked that the same is true for the other Higgs branch¬ 
ing ratios measured currently. For this reason, at the moment the Higgs branching ratios do 
not make important contributions to the global likelihood function of the CMSSM. However, 


10 As discussed later, the solid blue lines are the corresponding likelihood functions provided by the prospec¬ 
tive Higgs measurements at FCC-ee (TLEP). 
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Figure 17: The one-dimensional Ay 2 profile likelihood functions for the gluino mass nig 
(upper left panel), the generic first- and second-generation squark mass rrig (upper right 
panel), the lighter stop squark mass (lower left panel) and the lighter stau mass (lower right 
panel), as obtained using prospective FCC-ee (TLEP) precision electroweak measurements 
(solid red lines) and Higgs measurements (solid blue lines) J22 |/ with the same central values 
as the low-mass best-fit CMSSM point [l§j, neglecting the inevitable improvements in other 
constraints on the supersymmetric models. 


we see in the other panels of Fig. [18] that future measurements of the Higgs branching ra¬ 
tios at FCC-ee (TLEP) would have the potential to discriminate between different CMSSM 
parameter sets, so that much of the 68% and 95% CL regions in these panels are unshaded, 
since they lie more than three current standard deviations away from the prospective mea¬ 
surements. Specifically, several individual measurements at the central values predicted by 
the low-mass best-fit point in the CMSSM would each individually exclude regions at large 
values of mo and (particularly) m^. As in Fig. 15 for the electroweak precision observables, 


we see that prospective measurements of the observables studied are compatible with the 
low-mass best-fit values within one FCC-ee (TLEP) a only within fractions of the ‘Crimean’ 
68 % CL region. We also see that only narrow bands of the ‘Eurasian’ 95% CL regions would 
yield values of BR (H —>■ ZZ, 77) and BR (H —>■ gg) within one a of the low-mass best-fit 
prediction, and that the band for BR (H — » gg) does not overlap the others. 
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Figure 18: The present measurement of BR(H —» ZZ) (upper left panel), and prospective 
FCC-ee (TLEP) measurements J of BR(H —» ZZ) (upper right), BR(H —»■ 77 ,) (lower 
left) and BR(H —> gg) (lower right) are superposed on the (mo, mi/ 2 ) plane in the CMSSM 
shown previously in Fig. [1} The colours represent deviations from the present central value 
in units of the present LHC experimental error (upper left panel), and the deviations from 
the values at the low-mass best-fit CMSSM point (12] of the values at other points in the 


(mo,rai/ 2 ) plane in units of the estimated future FCC-ee (TLEP) experimental errors (22l 
(other panels). 
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Also as in the case of the electroweak precision observables discussed above, we have made 
a crude estimate of the impact of the prospective FCC-ee (TLEP) Higgs measurements 1 22 1 
on the global ;y 2 function for the CMSSM, again neglecting the inevitable improvements in 
flavour and dark matter observables, and setting aside the electroweak precision observables 
as well as the direct measurements of sparticle masses. As we see in Fig. 19, the high- 
precision Higgs measurements would, by themselves, again provide constraints on mo and 
ttt. 1 / 2 , which would be of comparable importance to those from the electroweak precision 
observables shown in the corresponding panel of Fig. [Tbl 



Figure 19: The Ay 2 = 2.30 (68% CL) and Ay 2 = 5.99 (95% CL) contours (red and blue, 
respectively) in the (m 0 i m i/ 2 ) plane for the CMSSM, assuming that Higgs measurements at 
FCC-ee (TLEP) ^22] have the same central values as at the current low-mass best-fit points, 
and neglecting inevitable improvements in other constraints on the supersymmetric models. 

We display as solid blue lines in Fig. [17] the corresponding one-dimensional projections of 
the contribution of prospective FCC-ee (TLEP) Higgs measurements to the CMSSM global 
X 2 function for m g (upper left panel), the generic squark mass (upper right panel), 
the lighter stop squark mass mq (lower left panel) and the lighter stau mass (lower 
right panel). We see that the Higgs estimates are comparable with the corresponding one¬ 
dimensional projections of the contribution of prospective FCC-ee (TLEP) electroweak pre¬ 
cision measurements, shown as solid lines in the various panels of Fig. HO The corresponding 
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95% CL mass ranges are estimated to be 


m-g e (1980,2500) GeV, 
rriq e (1780,2320) GeV, 
m fl e (370,510) GeV, 

m h e (890,1170) GeV, (12) 

whereas the nominal values at the best-ht point are 2280, 2080, 450 and 1020 GeV, respec¬ 
tively. The Higgs measurements clearly add another dimension to the tests of supersymmetric 
models at the loop level. 


4.4 Comparison between LHC and e + e Measurements 

The potential comparison between LHC and FCC-ee (TLEP) measurements in the best-fit 
low-mass CMSSM scenario can be seen in Fig. 20, where we overlay in the CMSSM (m 0 , m i/ 2 ) 
plane the potential direct measurements at the LHC presented earlier (pink and blue shad¬ 
ing) with indirect determinations at FCC-ee (TLEP) via EWPOs and Higgs measurements. 
A triple coincidence of direct sparticle mass measurements with indirect predictions from 
EWPOs and Higgs measurements would be truly impressive, a worthy successor to the suc¬ 
cessful predictions of the top and Higgs masses based on electroweak precision observables 
at LEP. 


4.5 Probes of Grand Unification 


The precision measurements of electroweak precision observables and the strong coupling 
ot s { m z) at LEP and the SLC opened a new chapter in probes of models of grand unifica¬ 
tion [541, making possible for the first time a clear discrimination between the predictions 
of supersymmetric and non-supersymmetric scenarios. As remarked in 122], it is clear that 
FCC-ee (TLEP) measurements could take this confrontation between experiment and differ¬ 
ent grand unified theories to a completely new level, through more accurate determinations 
of the SU(3), SU(2) and U(l) couplings by specifying the supersymmetric spectrum and 
hence TeV-scale threshold corrections to the running of the couplings. The combination 
of these measurements would enable powerful constraints to be placed on the GUT-scale 
particles in any specific GUT model. 

As an indication of this possibility, we consider the simplest supersymmetric SU(5) GUT, 
in which the GUT-scale particles comprise the heavy vector bosons V, the 24-plet Higgs 
bosons E and the coloured Higgs triplet bosons H c . By considering the three one-loop 
renormalization-group equations (RGEs) for the SU(3), SU(2) and U(l) couplings in this 
model, Murayama and Pierce (551 [^] derived the following one-loop relation between the 
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For a recent analysis, see 
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Figure 20: Upper panels: The 68 and 95% CL regions in the (m 0 , mi/ 2 ) plane of the CMSSM, 
overlaying potential direct measurements at LHC If with 300/fb (left panel) and 3000/fb 
(right panel) (pink and blue shading) with indirect determinations via electroweak precision 
and Higgs measurements at FCC-ee (TLEP) fli^ (red and blue solid lines). Lower pan¬ 
els: The 68% and 95% CL contours for the combination of the prospective constraints from 
LHC 14 with 300/fb (left panel) and 3000/fb (right panel) with the indirect determinations 
via electroweak precision and Higgs measurements at FCC-ee (TLEP) fH| / (red and blue solid 
lines). Also shown are the 68 and 95% CL regions found in a recent global fit to the CMSSM 
(red and blue dashed lines) [12]. 


low-energy values of the couplings, ai(m z ), the supersymmetric threshold and the mass of 
the Higgs triplet: 
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Another combination of the one-loop RGEs gives a similar relation for a combination of the 
masses of the heavy vector bosons V and the 24-plet Higgs bosons E: 
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These relations are subject to corrections from higher-order terms in the RGEs, etc., but may 
be used to estimate the uncertainties in the GUT-scale masses associated with uncertainties 
in the low-scale inputs. 


For example, inverting (13) we find 
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and inverting (14) we find 
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(16) 

It is estimated that at FCC-ee one could attain uncertainties A a^(m z ) ~ 10 -4 (correspond¬ 
ing to A a^ l {m z ) 10 2 ) and A sin 2 9w ~ 10 6 , with an input parametric uncertainty 
Att“^(m z ) ~ 5 x 10 -5 . Using 

1 sin 2 9w ^ 


a 2 (m z ) (m z ) ’ 


we find that 
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and infer that the dominant uncertainty in a 2 (m z ) is that due to A sin 2 9 W , giving us the 
estimate 

1 ' 10“ 4 . (19) 


Using then the relationship 
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it is evident that the dominant uncertainty in af 1 (mz) is correlated with that in a^mz): 
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( 21 ) 


It is apparent from the estimates (19) and (21) that the uncertainties due to the electroweak 
couplings in the GUT mass estimates (15) and (16) are much smaller than the uncertainties 
due to the strong coupling. 

Specifically, the precision measurements at FCC-ee (TLEP) should enable the mass of 
the colour-triplet to be estimated with an accuracy at the percent level, and similarly for 
the combination assuming that msusY can be determined with similar (or better) 

precision via direct or indirect measurements. Needless to say, this possibility of constraining 
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GUT-scale masses would apply within a specific GUT model, and the implications of the 
FCC-ee (TLEP) measurements would depend on the model. However, this analysis makes 
the point that high-precision measurements with FCC-ee (TLEP) could impose important 
constraints on GUT models, taking to the next level the insights provided previously by LEP 
measurements 1541. 


5 Prospects for the Discovery of Supersymmetry in 
Pessimistic Scenarios 

We now consider the prospects for discovering supersymmetry in ‘pessimistic’ high-mass 
CMSSM scenarios in which the HL-LHC does not discover supersymmetry, but provides 
only 95% CL lower limits on model parameters. 


5.1 Impact of LHC Searches 


To probe this case, we first make a crude estimate of the impact of such a negative result 
by including in the global ;y 2 functions for the CMSSM contributions based on the green 
lines in Fig. [lj which correspond to 95% CL exclusion by the LHC with 3000/fb of data, 
neglecting again the inevitable improvements in the measurements of other observables that 
would provide additional constraints on supersymmetry. The resulting 68% and 95% CL 
contours (red and blue, respectively) are shown in Fig. 21 The ‘Crimea’ region has now 


disappeared completely, and only the ‘Eurasia’ region remains. However, in the CMSSM, 
although this region is unified at the 95% CL, it is divided at the 68% CL into regions at 
lower and higher values of m 0 . 

shows the corresponding one-dimensional profile likelihood functions for nig (up- 


Fig. 


22 


per left panel) and nig (upper right panel). In the case of the gluino, we find a prospective 
95% CL lower limit nig > 3 TeV, and a lower limit rn, f] > 4 TeV for the squark mass. The 


limitations of our CMSSM sample 12 are such that we do not have any useful information 


about the likelihood functions for large masses, where they are expected to be quite flat. 
The lower left panel of Fig. [22] shows the corresponding one-dimensional profile likelihood 
function for rn^: the dip at m fl ~ 1000 GeV corresponds to the ‘cockscomb’ feature visible 


as an isolated 68% CL region with (mo,mi/ 2 ) ~ (1800,2400) GeV in Fig. 21, with the local 
peak at ~ 1300 GeV corresponding to the gap between the ‘cockscomb’ and the Eurasia 
region. As is also apparent in Fig. [13| within the CMSSM there is significant likelihood that 
nif x < 1500 GeV, so that f\ pair-production would be possible at CLIC, even if the LHC 
fails to discover supersymmetry. 
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Figure 21: The Ay 2 = 2.30 (68% CL) and Ax 2 = 5.99 (95% CL) contours (red and blue, 
respectively) in the (mo, mi/ 2 ) plane for the CMSSM, assuming that supersymmetry has not 
been discovered at the LHC with 3000/fb of luminosity, and neglecting inevitable improve¬ 
ments in other constraints on the supersymmetric models. 


5.2 Direct Sparticle Searches at a Higher-Energy Proton-Proton 
Collider 

We now turn to the potential of a future higher-energy hadron collider for discovering su¬ 
persymmetry within the CMSSM framework. To this end, we first analyze the nature of 
the CMSSM parameter space for large values of mo and mi/ 2 , taking into account the cold 
dark matter density constraint and the measurement of m/ ( ,, which are the only constraints 
capable of imposing upper limits on mo and mi/ 2 . Generally speaking, bringing the relic 
density down into the astrophysical range when these mass parameters are large requires 
some specific features in the sparticle spectrum such as near-degeneracy between the LSP, 
the NLSP and perhaps other supersymmetric particles, so as to suppress the relic LSP den¬ 
sity via coannihilation, or the existence of a massive Higgs boson that acts as an s-channel 
resonance and thereby suppresses the LSP density by enhancing LSP annihilation. 


One such possibility is the stau coannihilation strip 57 , which appears at low values of 


the ratio mo/mi/ 2 , adjacent to the stau LSP region. Its length depends on tan/5 and A 0 , 
extending as far as mi/ 2 ~ 1400 GeV for tan (5 = 40 and A 0 = 2.5m 0 |8|. A recent study 
has shown that all this strip may be explored by Run 2 of the LHC at 14 TeV |9|, as also 
discussed above. 

which appears at higher values of 


Another possibility is the focus-point strip 14,58 


the ratio mo/mi/ 2 , adjacent to the boundary of the region where one can find a consistent 
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Figure 22: The one-dimensional profile likelihood functions for rrig (upper left panel), rriq 
(upper right panel), (lower left panel) and rn^ (lower right panel) in the CMSSM, as¬ 
suming that supersymmetry has not been discovered at the LHC with 3000/fb of luminosity, 
and neglecting inevitable improvements in other constraints on the supersymmetric models. 


electroweak vacuum, and generally lies beyond the reach of the LHC searches discussed 
earlier. Along this strip, the Higgsino component of the neutralino LSP is enhanced, and 
its annihilations and coannihilations with heavier neutralinos and charginos are enhanced. 
Various studies have shown that the focus-point strip may extend to very large values of 
m 0 and with m 0 / m i /2 ~ 3 and A 0 < m 0 . The upper panels of Fig. 23 display a pair 


of focus-point strips for A 0 = 0 and tan/3 = 10 (left panel) and 52 (right panel), assuming 
m t = 173.2 GeV 12 . The regions of these planes where there is no consistent electroweak 


vacuum are coloured purple, the ochre regions at lower mo/mi /2 are excluded because of 
a charged LSP and/or a tachyon, and the green shaded regions are excluded by b —>■ sy 


12 


Larger (smaller) values of mo/rrii/2 are found for larger (smaller) values of mt. 
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13 


decay 

In both the upper panels of Fig. [23| we see a (dark blue) focus-point strip hugging 
the boundary of the region at mo/mi /2 > 3 where electroweak symmetry-breaking is not 
possible. In the right panel, we also see a rapid-annihilation funnel [4] projecting out of the 
stau coannihilation strip at low mo/mi /2 and extending to (mo,ffli/ 2 ) ~ (2500,1800) GeV. 
We also see in both panels that the contours of irq, calculated using FeynHiggs 2.10.0 [32] 
(shown as red dashed lines) are almost orthogonal to the focus-point strips. We note that the 
LHC measurement of even allowing for a 3 GeV uncertainty in the FeynHiggs 2.10.0 
calculation, excludes mp 2 > 4000 GeV. The solid black, blue, green and purple lines in each 
panel are particle exclusion reaches for Ifi T searches with LHC at 8 TeV, 300 and 3000/fb with 
LHC at 14 TeV, and 3000/fb with HE-LHC at 33 TeV, respectively. The focus-point strip 
extends beyond the reach of the LHC, even with 3000/fb at 14 TeV in the centre of mass, 
and even beyond the reach of the HE-LHC with 3000/fb at 33 TeV. However, the portion 
allowed by the Higgs mass constraint lies comfortably within the reach of the FCC-hh with 
3000/fb at 100 TeV, as discussed below. 

at larger values of A 0 /m 0 >2.2 


As shown by two examples in the lower panels of Fig. 


23 


there are wedges at larger mo/mi /2 that are excluded because the lighter stop squark is the 
LSP. Hugging the boundaries of these wedges there are narrow stop coannihilation strips 115, 


59 


where 8m q — m x is small. 


The opening angles of the stop LSP wedges have little 
dependence on tan/3, and both the planes we show have tan/3 = 20. On the other hand, 
the opening angles of the stop LSP wedges increase with A 0 /m 0 , with the result that the 
wedge at intermediate mo/mi /2 where the LSP is the lightest neutralino is closed off for 
A 0 /m 0 > 5.5. 


We display in the lower panels of Fig. 23 the cases tan/3 = 2.3 (left panel) and 3 (right 
panel). As we discuss later, the stop coannihilation strips extend far beyond the ranges of mo 


and mi /2 that we display in these two panels of Fig. 23. These panels also display contours of 
mf L calculated using FeynHiggs 2.10.0. We note that, in contrast to the focus-point cases 
displayed in Fig. [23j the irq, contours cross the dark matter strip at a much smaller angle. 
As a consequence, the allowed ranges of mo and mi/ 2 are much larger than in the focus- 
point case, after allowing for the 3 GeV uncertainty estimated within FeynHiggs 2.10.0 
for a given input SLHA file 1601. Additionally, we find that the value irq, calculated along 
the stop coannihilation strip is very sensitive to the codes used to evolve the CMSSM input 
parameters down to low energies and calculate the spectra used as SLHA inputs in the 
FeynHiggs 2.10.0 calculation of ?rq t . This introduces an additional uncertainty of several 
GeV in the value of /rq, corresponding to any given set of CMSSM input parameters. For 
this reason, no portions of these stop coannihilation strips can be excluded on the basis of 
the LHC measurement of /rq, 


14 


The lighter green line in the upper right panel of Fig. 23 delineates the region allowed by B s —» //+// 


We assume in calculating these parameter planes that a s (mz ) = 0.1184 and sin" f)\\- = 0.2325. With 
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Figure 23: Upper panels: the (mo,mi/ 2 ) planes for A 0 = 0 and tan/3 = 10 (left panel) 
and 40 (right panel), displaying focus-point strips. Lower panels: the (mo,mi/ 2 ) planes for 
tan/3 = 20 and Aq/tuq = 2.3 (left panel), and A 0 /m 0 = 3.0 (right panel). In each panel, 
the ochre regions are excluded because of a charged LSP and/or a tachyon, and the green 
regions are excluded by b —>• sy decay. There is no consistent electroweak vacuum in the 
purple regions in the upper panels. In the dark blue strips the relic LSP density lies within 
the range allowed by cosmology, and the dashed red lines are contours of nth as calculated 
using FeynHiggs 2.10.0. The solid black, blue, green and purple lines in each panel are 
particle exclusion reaches fortyr searches with LHC at 8 Te V, 300 and 3000/fb with LHC at 
If TeV, and 3000/fb with HE-LHC at 33 TeV, respectively. 


these choices, the SU(3), SU(2) and U(l) couplings are in general not quite equal at the GUT scale, but the 
deviations are generally compatible with zero when 3late uncertainties in the inputs are taken into account, 
and any apparent non-universality could be cancelled by GUT threshold corrections and/or produced by 
higher-dimensional operator contributions. The comparable planes shown in 15 assumed sin 2 9w = 0.2325 


and strict coupling-constant unification, and consequently had a range of values of a s (mz ). This is the 
reason for the differences between these planes and comparable planes in 
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As in the upper panels of Fig. [23| the solid black, blue, green and purple lines in the lower 
panels are particle exclusion reaches for Ifir searches with LHC at 8 TeV, 300 and 3000/fb 
with LHC at 14 TeV, and 3000/fb with HE-LHC at 33 TeV, respectively, now for CMSSM 
scenarios with stop coannihilation strips. We see that the LHC sensitivity contours in the 
lower panels of Fig. [23] include only portions of these stop coannihilation strips extending to 
my 2 ~ 500 (700) GeV for tan /3 = 20 and A 0 = 2.3 (3) mo. The HE-LHC sensitivity contour 
with 3000/fb at 33 TeV extends to ~ 2400(3500) GeV along the stop coannihilation 
strips for tan/3 = 20 and A 0 = 2.3 (3) mo, but large fractions of these strips lie beyond its 
reach. 

Fig. pijFj displays the profiles of the focus-point strips in Fig. 23 (upper panels) and 


of the stop coannihilation strips in Fig. 23 (lower panels), along their full lengths. Both 
pairs of profiles exhibit the values of m/ ( , calculated using SLHA files obtained using SSARD 
as inputs to FeynHiggs 2.10.0 (near-horizontal solid green fines), including uncertainty 
estimates of ±3 GeV (near-horizontal dashed green lines). As already noted, only portions 
of the focus-point strips are compatible with the LHC measurement of nit, (yellow bands) 
within these uncertainties, whereas in the cases of the stop coannihilation strips there are 
significant additional uncertainties associated with the RGE running, and all portions of the 
strips are compatible with m/ t . In the cases of the stop coannihilation strips in the lower 
panels of Fig. 24, we also display as blue lines the mass difference 5m = my — rn x along the 
strips 16 . In the examples shown, this mass difference is generally < mw + mb, so that the 


branching ratio for two-body t\ —? \ + c decay usually dominates over that for four-body 
t\ —yx + W + b + v decay. However, this is not always the case, as illustrated by examples 


in 115] and by Fig. 25 for the stop coannihilation strip with A 0 /m 0 = 3.0 and tan/3 = 20. 
The branching ratio for t\ —» x + W + b + v decay may dominate when uy — m x > mw + mb, 
as seen in the lower right panel of Fig. [24| Thus, a complete search for supersymmetry at 
FCC-hh should include searches for both the t\ —> x + c and i\ X + H 7 A b + v decay 
signatures. 

The (near-)vertical lines in Fig. 24 mark our estimates of the sensitivities of the LHC 
(black - 8 TeV, blue - 300/fb at 14 TeV, green - 3000/fb at 14 TeV), 3000/fb at HE-LHC 
(purple) and 3000/fb at FCC-hh (red) along the stop coannihilation strips. The solid lines 
represent the extrapolated reaches of the generic jets + 1£t searches, and the dashed lines 
in the lower panels represent the extrapolated reaches of dedicated searches for t\ —> c + x 
decays, which lose some sensitivity as 5m increases because of the increase in the t\ —y 
X + W + b + v decay branching ratio. We see that the FCC-hh would be sensitive to the 
full extents of the focus-point strips (upper panels) and of the stop coannihilation strip for 
A 0 = 2.3 mo (lower left panel), but not all the stop coannihilation strip for A 0 = 3.0 mo 


24 


15 The upper panels of Fig. 

16 This is generally larger than in the cases with strict gauge coupling unification assumed in 15 


in the published version are corrected in an erratum 


corollary that the strips extend to larger m 1/2 than shown there. 


with the 
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Figure 24: Upper panels: the solid blue lines are the profiles in the (m 0 , plane of the 

focus-point strips for A$ = 0 and tan/3 = 10 (left panel), and A 0 = 0 and tan /3 = 52 (right 
panel). Lower panels: the solid blue lines are the profiles in the (mi/ 2 , Sm = —m x ) plane 
of the stop coannihilation strips for Aq/uiq = 2.3 and tan j3 — 20 (left panel), and Aq/uiq = 
3.0 and tan (3 = 20 (right panel). The near-vertical black, blue, green, purple and red lines in 
each panel are particle exclusion reaches for particle searches with LHC at 8 Te V, 300 and 
3000/fb with LHC at If TeV, 3000/fb with HE-LHC at 33 TeV and 3000/fb with FCC-hh at 
100 TeV, respectively. The solid lines are for genericIpT searches, and (in the lower panels) 
the dashed lines are for dedicated stop searches. The solid (dashed) near-horizontal green 
lines are central values (probable ranges) of mn calculated using FeynHiggs 2.10.0, and the 
yellow band represents the experimental value of 
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tan/5 = 20, Cl x h 2 = 0.12,p > 0 



Figure 25: The branching ratio for t\ —» y + c along the stop coannihilation strips for 
tan/5 = 20 and A 0 /mo = 2.3 (solid black line) and Aq/tuq = 3.0 (dashed blue line). In the 
latter case the branching ratio drops to a minimum <0.1 when — m x > mw + fnb, as 
seen in the lower right panel of Fig. \2f\ 


(lower right panel): this is true in general for Aq/tuq > 2.5. 


5.3 Impacts of Electroweak and Higgs Precision Observables 


We have also studied the possible impact of EWPOs and Higgs precision measurements 
along the focus-point and stop coannihilation strips discussed in the previous Subsection. 
We find that the contributions to the global y 2 function of the present EWPOs and Higgs 
measurements do not vary strongly along the strips, so do not discuss them further. Instead, 
we focus on the potential impacts of FCC-ee measurements, choosing benchmark points 
on these strips. These benchmark points are chosen to have values of run, as calculated 
with FeynHiggs 2.10.0, that are highly compatible with the central experimental value 
m h ~ 125 GeV0 

Fig. [26| shows the estimated contributions of the EWPOs and Higgs observables measured 
at FCC-ee (TLEP) (red and blue lines, respectively) to the global y 2 functions along the 
focus-point and stop coannihilation strips, which are plotted using mo as the horizontal axis. 
The diagonal dashed black lines show the corresponding values of using the scale shown 


17 We recall that calculations of ///with FeynHiggs 2.10.0 have uncertainties ~ ±3 GeV so that, as 
substantial fractions of the focus-point and stop coannihilation strips are compatible with 


discussed in 15 


the mh constraint. 
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on the right-hand vertical axis. In each case, we have assumed measurements with FCC- 
ee (TLEP) uncertainties and central values coinciding with those calculated using CMSSM 
model parameters at the benchmark points shown as the black spots. In each case, we see 
that the FCC-ee (TLEP) measurements would be capable of specifying with an accuracy 
that is greatest for the stop coannihilation strip with A 0 /mo = 3.0 and tan/5 = 20 (lower 
right panel) and least for the focus-point strip with A 0 = 0 and tan (3 = 10 (upper left panel). 
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Figure 26: Upper panels: the estimated contributions of the EWPOs and Higgs observables 
measured at FCC-ee (TLEP) to the global y 2 function (solid red and blue lines, respectively) 
assuming the parameters of benchmark points along the focus-point strips for A 0 = 0 and 
tan/5 = 10 (left panel), and Aq = 0 and tan /5 = 52 (right panel). Lower panels: the same for 
the stop coannihilation strips for A 0 /m 0 = 2.3 and tan /5 = 20 (left panel), and A 0 /m 0 = 3.0 
and tan/5 = 20 (right panel). The diagonal black dashed lines show the values of m. 1/2 along 
the corresponding strips (right-hand vertical axes), and the black spots show the parameters 
of the corresponding benchmark points. 

All of these benchmark points he within reach of generic Ifr searches at FCC-hh, as 
seen in Fig. [24} (Indeed, the fixed-point benchmark points he with the reach of HE-LHC.) 
Therefore, for each of these benchmark points it would be possible to make a comparison be¬ 
tween the direct determination of the CMSSM parameters with those inferred indirectly from 
FCC-ee measurements, much as we discussed earlier for the case of the low-mass CMSSM 
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benchmark point, the LHC and FCC-ee. In all cases, FCC-ee would make possible tests 
of supersymmetry at the loop level, even in pessimistic scenarios where the LHC does not 
discover supersymmetry. 


6 Summary 

We have explored in this paper the interplay between direct and indirect searches for su¬ 
persymmetry in future runs of the LHC and at proposed future colliders. This is clearly a 
very broad topic, so we have restricted our attention to the CMSSM. A recent global fit to 
the CMSSM has found two favoured regions of its parameter space: a low-mass ‘Crimea’ 
region and a high-mass ‘Eurasia’ region. In the optimistic low-mass case, extrapolating the 
sensitivities of supersymmetry searches with LHC Run 1, we have found that future runs 
of the LHC with 300 or 3000/fb of data at 14 TeV should be able to discover gluinos and 
squarks if Nature is described by the CMSSM in the Crimea region. Moreover, the LHC 
experiments should be able to measure the gluino and squark masses with high accuracy and 
hence also the soft supersymmetry-breaking parameters of the CMSSM. 

In this optimistic scenario, electroweakly-interacting sparticles could be discovered at 
the CLIC e + e _ collider and their masses measured very accurately, providing an important 
cross-check on the consistency of the CMSSM with its universal input soft supersymmetry¬ 
breaking parameters. On the other hand, an e + e _ collider with centre-of-mass energy limited 
to 0.5 TeV would not produce supersymmetric particles even in this optimistic scenario, 
though a collider with a centre-of-mas energy of 1 TeV could detect the lighter stau slepton. 

However, an e + e _ collider capable of high-luminosity running on the Z peak as well 
as producing large numbers of Higgs bosons, such as FCC-ee (TLEP), would provide two 
sets of high-precision measurements that could be used to constrain supersymmetric loop 
corrections and hence, in conjunction with the LHC measurements, check supersymmet¬ 
ric predictions at the quantum level in two independent ways. The combination of direct 
and indirect measurements possible in this optimistic scenario would test the CMSSM in a 
way reminiscent of the use of precision measurements from LEP and elsewhere to predict 
successfully the masses of the top quark and the Higgs boson, albeit in time-reversed order. 

In the pessimistic scenario in which sparticles are too heavy to be produced at the LHC, 
the CLIC e + e _ collider might still be able to discover the lighter stau slepton. We have also 
assessed the capability of a higher-energy proton-proton collider such as FCC-hh to discover 
supersymmetry, and the ability of high-precision FCC-ee (TLEP) measurements to provide 
any indirect evidence. For this part of our analysis, we concentrate on the narrow strips 
in the CMSSM parameter space that extend to large sparticle masses, namely the stop- 
coannihilation strip and the focus-point strip. We find that a 33-TeV collider such as HE- 
LHC would be able to discover supersymmetry via I£t searches along some fractions of these 
strips, and that a 100-TeV collider such as FCC-hh would be able to discover supersymmetry 
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along most of the extents of the strips examined. By studying specific benchmark points 
along these strips, we have shown that FCC-ee measurements could in principle determine 
indirectly CMSSM model parameters also in the pessimistic scenario, and the combination 
with FCC-hh measurements could test supersymmetry at the loop level. 

In both the optimistic and pessimistic scenarios, we find that high-precision measure¬ 
ments with a high-luminosity e + e _ collider can play roles that are complementary to direct 
particle searches with proton-proton colliders, and could enable supersymmetry to be tested 
at the quantum (loop) level). Run 2 of the LHC will provide us with some valuable pointers 
indicating which of these scenarios may be realised in Nature. 


Acknowledgements 

The work of K.J.dV., O.B. and J.E. is supported in part by the London Centre for Terauni- 
verse Studies (LCTS), using funding from the European Research Council via the Advanced 
Investigator Grant 267352. The work of J.E. is also supported in part by STFC (UK) via 
the research grants ST/J002798/1 and ST/L000326/1. The work of K.A.O. is supported in 
part by DOE grant de-sc0011842 at the University of Minnesota. 


References 

[1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716 , 1 (2012) |arXiv:1207.7214 
[hep-ex]]; S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716 , 30 (2012) 
|arXiv: 1207.7235 [hep-ex]]. 

[2] G. Aad et al. [ATLAS Collaboration], JHEP 1409 (2014) 176 |arXiv:1405.7875 
[hep-ex]]; full ATLAS Run 1 results can be found at 

https://twiki.cern.ch/twiki/bin/view/AtlasPublic/SupersymmetryPublicResults. 

[3] S. Chatrchyan et al. [CMS Collaboration], JHEP 1406 (2014) 055 

|arXiv: 1402.4770 [hep-ex]]; full CMS Run 1 results can be found at 

https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsSUS. 

[4] M. Drees and M. M. Nojiri, Phys. Rev. D 47 (1993) 376 |arXiv:hep-ph/9207234|; H. Baer 
and M. Brhlik, Phys. Rev. D 53 (1996) 597 |arXiv:hep-ph/9508321|; Phys. Rev. D 
57 (1998) 567 |arXiv:hep-ph/9706509|; H. Baer, M. Brhlik, M. A. Diaz, J. Ferrandis, 

P. Mercadante, P. Quintana and X. Tata, Phys. Rev. D 63 (2001) 015007 |arXiv:hep- 
ph/0005027|; J. R. Ellis, T. Falk, G. Ganis, K. A. Olive and M. Srednicki, Phys. Lett. 

B 510 , 236 (2001) |hep-ph/0102098]. 

[5] G. L. Kane, C. F. Kolda, L. Roszkowski and J. D. Wells, Phys. Rev. D 49 (1994) 6173 
|arXiv:hep-ph/9312272|; J. R. Ellis, T. Falk, K. A. Olive and M. Schmitt, Phys. Lett. 


41 


B 388 (1996) 97 |arXiv:hep-ph/9607292|; Phys. Lett. B 413 (1997) 355 |arXiv:hep- 
ph/9705444|; J. R. Ellis, T. Falk, G. Ganis, K. A. Olive and M. Schmitt, Phys. Rev. 
D 58 (1998) 095002 |arXiv:hep-ph/9801445|; V. D. Barger and C. Kao, Phys. Rev. D 
57 (1998) 3131 |arXiv:hep-ph/9704403|; J. R. Ellis, T. Falk, G. Ganis and K. A. Olive, 
Phys. Rev. D 62 (2000) 075010 |arXiv:hep-ph/0004169|; L. Roszkowski, R. Ruiz de Aus- 
tri and T. Nihei, JHEP 0108 (2001) 024 |arXiv:hep-ph/0106334|; A. Djouadi, M. Drees 
and J. L. Kneur, JHEP 0108 (2001) 055 |arXiv:hep-ph/0107316|; U. Chattopadhyay, 
A. Corsetti and P. Nath, Phys. Rev. D 66 ( 2002 ) 035003 |arXiv:hep-ph/0201001|; 

J. R. Ellis, K. A. Olive and Y. Santoso, New Jour. Phys. 4 ( 2002 ) 32 |arXiv:hep- 
ph/ 0202110 ]. 

[ 6 ] H. Goldberg, Phys. Rev. Lett. 50 (1983) 1419; J. Ellis, J. Hagelin, D. Nanopoulos, 

K. Olive and M. Srednicki, Nucl. Phys. B 238 (1984) 453. 

[7] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Lett. B 565 (2003) 176 

|arXiv:hep-ph/0303043|; H. Baer and C. Balazs, JCAP 0305, 006 (2003) |arXiv:hep- 

ph/0303114||; A. B. Lahanas and D. V. Nanopoulos, Phys. Lett. B 568, 55 (2003) 
|arXiv:hep-ph/0303130|; U. Chattopadhyay, A. Corsetti and P. Nath, Phys. Rev. D 68 , 
035005 (2003) |arXiv:hep-ph/030320l]; C. Munoz, Int. J. Mod. Phys. A 19, 3093 (2004) 
|arXiv:hep-ph/0309346|; R. Arnowitt, B. Dutta and B. Hu, arXiv:hep-ph/0310103; 
J. Ellis and K. A. Olive, in Particle dark matter, ed. G. Bertone (Cambridge University 
Press, 2010 ) ppl42-163 |arXiv:1001.3651 [astro-ph.CO]]. 

[ 8 ] J. Ellis and K. A. Olive, Eur. Phys. J. C 72, 2005 ( 2012 ) |arXiv:1202.3262 [hep-ph]]. 

[9] O. Buchmueller, M. J. Dolan, J. Ellis, T. Hahn, S. Heinemeyer, W. Hollik, J. Marrouche 
and K. A. Olive et al ., Eur. Phys. J. C 74 (2014) 3, 2809 |arXiv:1312.5233 [hep-ph]]. 

[10] O. Buchmueller, et al, Eur. Phys. J. C 72 ( 2012 ) 2020 |arXiv:1112.3564 [hep-ph]]. 

[ 11 ] O. Buchmueller, et al, Eur. Phys. J. C 72 ( 2012 ) 2243 |arXiv:1207.7315 [hep-ph]]. 

[ 12 ] O. Buchmueller, R. Cavanaugh, A. De Roeck, M. J. Dolan, J. R. Ellis, H. Flacher, 
S. Heinemeyer and G. Isidori et al, Eur. Phys. J. C 74 (2014) 6 , 2922 |arXiv:1312.5250 
[hep-ph]]. 

[13] H. Baer, V. Barger and A. Mustafayev, Phys. Rev. D 85, 075010 ( 2012 ) |arXiv:1112.3017 
[hep-ph]]; T. Li, J. A. Maxin, D. V. Nanopoulos and J. W. Walker, Phys. Lett. B 710 
( 2012 ) 207 |arXiv:1112.3024 [hep-ph]]; S. Heinemeyer, O. Stal and G. Weiglein, Phys. 
Lett. B 710, 201 ( 2012 ) jarXiv: 1112.3026 [hep-ph]]; A. Arbey, M. Battaglia, A. Djouadi, 
F. Mahmoudi and J. Quevillon, Phys. Lett. B 708 ( 2012 ) 162 |arXiv:1112.3028 [hep- 
ph]]; P. Draper, P. Meade, M. Reece and D. Shih, Phys. Rev. D 85, 095007 ( 2012 ) 


42 



|arXiv:1112.3068 [hep-ph]]; S. Akula, B. Altunkaynak, D. Feldman, P. Nath and 
G. Peim, Phys. Rev. D 85 (2012) 075001 |arXiv:1112.3645 [hep-ph]]; M. Kadastik, 
K. Kannike, A. Racioppi and M. Raidal, JHEP 1205 (2012) 061 |arXiv:1112.3647 
[hep-ph]]; C. Strege, G. Bertone, D. G. Cerdeno, M. Fornasa, R. Ruiz de Austri and 
R. Trotta, JCAP 1203, 030 (2012) |arXiv:1112.4192 [hep-ph]]; J. Cao, Z. Heng, D. Li 
and J. M. Yang, Phys. Lett. B 710 (2012) 665 jarXiv: 1112.4391 [hep-ph]]; L. Apari- 
cio, D. G. Cerdeno and L. E. Ibanez, JHEP 1204, 126 (2012) |arXiv:1202.0822 [hep- 
ph]]; H. Baer, V. Barger and A. Mustafayev, JHEP 1205 (2012) 091 |arXiv:1202.4038 
[hep-ph]]; P. Bechtle, T. Bringmann, K. Desch, H. Dreiner, M. Hamer, C. Hensel, 
M. Kramer and N. Nguyen et al, JHEP 1206, 098 (2012) jarXiv: 1204.4199 [hep-ph]]; 
C. Balazs, A. Buckley, D. Carter, B. Farmer and M. White, arXiv:1205.1568 [hep- 
ph]; D. Ghosh, M. Guchait, S. Raychaudhuri and D. Sengupta, Phys. Rev. D 86, 
055007 (2012) |arXiv: 1205.2283 [hep-ph]]; A. Fowlie, M. Kazana, K. Kowalska, S. Mu¬ 
nir, L. Roszkowski, E. M. Sessolo, S. Trojanowski and Y. -L. S. Tsai, Phys. Rev. D 86, 
075010 (2012) |arXiv: 1206.0264 [hep-ph]]; K. Kowalska et al. [BayesFITS Group Collab¬ 
oration], Phys. Rev. D 87, no. 11, 115010 (2013) |arXiv: 1211.1693 [hep-ph]]; C. Strege, 
G. Bertone, F. Feroz, M. Fornasa, R. Ruiz de Austri and R. Trotta, JCAP 1304, 013 
(2013) |arXiv: 1212.2636 [hep-ph]]; M. E. Cabrera, J. A. Casas and R. R. de Austri, 
JHEP 1307 (2013) 182 [arXiv:1212.4821 [hep-ph]]; T. Cohen and J. G. Wacker, JHEP 
1309 (2013) 061 |arXiv:1305.2914 [hep-ph]]; S. Henrot-Versill, Rm. Lafaye, T. Plehn, 
M. Rauch, D. Zerwas, S. p. Plaszczynski, B. Rouill d’Orfeuil and M. Spinelli, Phys. Rev. 
D 89, 055017 (2014) |arXiv: 1309.6958 [hep-ph]]; P. Bechtle, K. Desch, H. K. Dreiner, 
M. Hamer, M. Krer, B. O’Leary, W. Porod and X. Prudent et al, arXiv: 1310.3045 
[hep-ph]; L. Roszkowski, E. M. Sessolo and A. J. Williams, arXiv: 1405.4289 [hep-ph], 

[14] J. L. Feng, K. T. Matchev and D. Sanford, Phys. Rev. D 85, 075007 (2012) 
|arXiv:1112.3021 [hep-ph]]; P. Draper, J. Feng, P. Kant, S. Profumo and D. Sanford, 
Phys. Rev. D 88, 015025 (2013) [arXiv: 1304.1159 [hep-ph]]. 

[15] J. Ellis, K. A. Olive and J. Zheng, Eur. Phys. J. C 74 (2014) 2947 |arXiv:1404.5571 
[hep-ph]]. 

[16] M. Battaglia, J. J. Blaising, J. S. Marshall, S. Poss, A. Sailer, M. Thomson and E. van 
der Kraaij, JHEP 1309 (2013) 001 jarXiv: 1304.2825 [hep-ex]]. 

[17] P. Lebrun, L. Linssen, A. Lucaci-Timoce, D. Schulte, F. Simon, S. Stapnes, N. Toge 
and H. Weerts et al, The CLIC Programme: Towards a Staged e + e~ Linear Collider 
Exploring the Terascale : CLIC Conceptual Design Report , arXiv: 1209.2543 [physics.ins- 
det]. 


43 


[18] S. Schael et al. [ALEPH, DELPHI, L3, OPAL and SLD Collaborations, LEP Elec- 
troweak Working Group and SLD Electroweak and Heavy Flavour Groups], Phys. Rept. 
427 (2006) 257 |hep-ex/0509008|. 

[19] ATLAS Collaboration, https://twiki.cern.ch/twiki/bin/ 
view/AtlasPublic/HiggsPublicResults. 

[20] S. Chatrchyan et al. [CMS Collaboration], JHEP 1306 (2013) 081 |arXiv:1303.4571 
[hep-ex]]. 

[21] CDF and DO Collaborations, http://tevnphwg.fnal.gov. 

[22] M. Bicer et al. [TLEP Design Study Working Group Collaboration], JHEP 1401 (2014) 
164 |arXiv: 1308.6176 [hep-ex]]. 

[23] https://espace2013.cern.ch/fcc/Pages/default.aspx. 

[24] For some studies with similar motivations, see A. Fowlie and M. Raidal, Eur. Phys. J. 
C 74 (2014) 2948 |arXiv: 1402.5419 [hep-ph]]; M. Low and L. T. Wang, JHEP 1408 
(2014) 161 |arXiv: 1404.0682 [hep-ph]]; B. S. Acharya, K. Bozek, C. Pongkitivanichkul 
and K. Sakurai, JHEP 1502 (2015) 181 |arXiv:1410.1532 [hep-ph]]; S. Gori, S. Jung, 

L. T. Wang and J. D. Wells, JHEP 1412 (2014) 108 [arXiv:1410.6287 [hep-ph]]; J. Bra- 
mante, P. J. Fox, A. Martin, B. Ostdiek, T. Plehn, T. Schell and M. Takeuchi, Phys. 
Rev. D 91 (2015) 5, 054015 |arXiv:1412.4789 [hep-ph]]. For some other studies of super- 
symmetry at a 100-TeV pp collider, see T. Cohen, T. Golling, M. Hance, A. Henrichs, 
K. Howe, J. Loyal, S. Padhi and J. G. Wacker, JHEP 1404 (2014) 117 |arXiv:1311.6480 
[hep-ph]]; T. Cohen, R. T. D’Agnolo, M. Hance, H. K. Lou and J. G. Wacker, JHEP 
1411 (2014) 021 |arXiv: 1406.4512 [hep-ph]]; C. Borschensky, M. Krmer, A. Kulesza, 

M. Mangano, S. Padhi, T. Plehn and X. Portell, Eur. Phys. J. C 74 (2014) 12, 3174 
|arXiv: 1407.5066 [hep-ph]]; S. A. R. Ellis, G. L. Kane and B. Zheng, arXiv:1408.1961 
[hep-ph]; H. Beauchesne, K. Earl and T. Gregoire, arXiv: 1503.03099 [hep-ph], 

[25] H. Baer, A. Mustafayev, S. Profumo, A. Belyaev and X. Tata, Phys. Rev. D 71, 095008 
(2005) [arXiv:hep-ph/0412059] and JHEP 0507 (2005) 065, hep-ph/0504001. 

[26] J. R. Ellis, K. A. Olive and P. Sandick, Phys. Rev. D 78, 075012 (2008) |arXiv:0805.2343 
[hep-ph]]; J. Ellis, F. Luo, K. A. Olive and P. Sandick, Eur. Phys. J. C 73, 2403 (2013) 
|arXiv: 1212.4476 [hep-ph]]. 

[27] J. Ellis, K. Olive and Y. Santoso, Phys. Lett. B 539, 107 (2002) |arXiv:hep-ph/0204192|; 
J. R. Ellis, T. Falk, K. A. Olive and Y. Santoso, Nucl. Phys. B 652, 259 (2003) 
|arXiv:hep-ph/0210205|. 


44 



[28] O. Buchmueller, R. Cavanaugh, M. Citron, A. De Roeck, M. J. Dolan, J. R. Ellis, 
H. Flaecher and S. Heinemeyer et al, Eur. Phys. J. C 74 (2014) 3212 |arXiv:1408.4060 
[hep-ph]]. 

[29] K. J. de Vries, E. A. Bagnaschi, O. Buchmueller, R. Cavanaugh, M. Citron, A. De 
Roeck, M. J. Dolan and J. R. Ellis et al, arXiv: 1504.03260 [hep-ph], 

[30] G. Aad et al. [CMS Collaboration], arXiv: 1503.07589 [hep-ex]. 

[31] G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich and G. Weiglein, Eur. Phys. J. 
C 28 (2003) 133 |arXiv:hep-ph/0212020|; S. Heinemeyer, W. Hollik and G. Wei¬ 
glein, Eur. Phys. J. C 9 (1999) 343 |arXiv:hep-ph/9812472|; S. Heinemeyer, W. Hol- 
lik and G. Weiglein, Comput. Phys. Commun. 124 (2000) 76 |arXiv:hep-ph/9812320|; 
M. Frank et al. , JHEP 0702 (2007) 047 |arXiv:hep-ph/0611326|; T. Hahn, S. Heine¬ 
meyer, W. Hollik, H. Rzehak and G. Weiglein, Comput. Phys. Commun. 180 (2009) 
1426. see http://www.feynhiggs.de . 

[32] T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak and G. Weiglein, |arXiv:1312.4937 [hep- 
ph]]. 

[33] [The Muon g-2 Collaboration], Phys. Rev. Lett. 92 (2004) 161802, hep-ex/0401008; 
G. Bennett et al. [The Muon g-2 Collaboration], Phys. Rev. D 73 (2006) 072003 
|arXiv:hep-ex/0602035|. 

[34] S. Chen et al. [CLEO Collaboration], Phys. Rev. Lett. 87 (2001) 251807 |arXiv:hep- 
ex/0108032|; P. Koppenburg et al. [Belle Collaboration], Phys. Rev. Lett. 93 (2004) 
061803 |arXiv:hep-ex/0403004|; B. Aubert et al. [BaBar Collaboration], arXiv:hep- 
ex/0207076; E. Barberio et al. [Heavy Flavor Averaging Group (HFAG)], arXiv:hep- 
ex/0603003. 

[35] V. M. Abazov et al. [DO Collaboration], Phys. Lett. B 693 (2010) 539 |arXiv:1006.3469 
[hep-ex]]; T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 107 (2011) 191801 
[Publisher-note 107 (2011) 239903] |arXiv:1107.2304 [hep-ex]]; G. Aad et al. [ATLAS 
Collaboration], Phys. Lett. B 713 (2012) 387 |arXiv: 1204.0735 [hep-ex]]. 

[36] R.Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. Ill (2013) 101805 

|arXiv: 1307.5024 [hep-ex]]. 

[37] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. Ill (2013) 101804 
|arXiv: 1307.5025 [hep-ex]]. 

[38] V. Khachatryan et al. [CMS and LHCb Collaborations], DOI: 10.1038/naturel4474, 
arXiv: 1411.4413 [hep-ex]. 


45 




[39] D. S. Akerib et al [LUX Collaboration], Phys. Rev. Lett. 112 , 091303 (2014) 
|arXiv: 1310.8214 [astro-ph.CO]]. 

[40] P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 571 (2014) A16 
|arXiv: 1303.5076 [astro-ph.CO]]. Using the Planck 2015 results: P. A. R. Ade et al. 
[Planck Collaboration], arXiv: 1502.01589 [astro-ph.CO] would have no visible effect on 
our analysis. 

[41] Y. Gershtein, M. Luty, M. Narain, L.-T. Wang, D. Whiteson, K. Agashe, L. Apanasevich 
and G. Artoni et al., Working Group Report: New Particles, Forces, and Dimensions, 
arXiv:1311.0299 [hep-ex], 

[42] CMS Collaboration, arXiv:1307.7135 [hep-ex]. 

[43] S. Malik, C. McCabe, H. Araujo, A. Belyaev, C. Boehm, J. Brooke, O. Buchmueller 
and G. Davies et al, arXiv: 1409.4075 [hep-ex], 

[44] O. Buchmueller, M. J. Dolan, S. A. Malik and C. McCabe, JHEP 1501 (2015) 037 
|arXiv: 1407.8257 [hep-ph]]; 

[45] G. Salam and A. Weiler, http://collider-reach.web.cern.ch/collider-reach/. 

[46] ATLAS Collaboration, https://cdsweb.cern.ch/record/1472518/files/ 
ATL-PHYS-PUB-2012-001.pdf. 

[47] T. Sjostrand, S. Mrenna and P. Z. Skands, Comput. Phys. Commun. 178 (2008) 852 
|arXiv:0710.3820 [hep-ph]]. 

[48] N. Desai and P. Z. Skands, Eur. Phys. J. C 72, 2238 (2012) |arXiv:1109.5852 [hep-ph]]. 

[49] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63 (2009) 189 
|arXiv:0901.0002 [hep-ph]]. 

[50] R. D. Ball, V. Bertone, S. Carrazza, C. S. Deans, L. Del Debbio, S. Forte, A. Guffanti 
and N. P. Hartland et al, Nucl. Phys. B 867 (2013) 244 |arXiv: 1207.1303 [hep-ph]], 
https://nnpdf.hepforge.org. 

[51] C. G. Lester and D. J. Summers, Phys. Lett. B 463 (1999) 99 |hep-ph/9906349|; A. Barr, 
C. Lester and P. Stephens, J. Phys. G 29 (2003) 2343 |hep-ph/0304226|. 

[52] M. Baak et al. [Glitter Group], Eur. Phys. J. C 74 (2014) 3046 |arXiv: 1407.3792 [hep- 
ph]]. 


46 



[53] A. Blondel, A. Chao, W. Chou, J. Gao, D. Schulte and K. Yokoya, Report of the 
ICFA Beam Dynamics Workshop ’Accelerators for a Higgs Factory: Linear vs. Circular’ 
(HF2012), arXiv:1302.3318 [physics, acc-ph], 

[54] J. Ellis, S. Kelley and D.V. Nanopoulos, Phys. Lett. B 249 (1990) 441; Phys. Lett. B 
260 (1991) 131; U. Amaldi, W. de Boer and H. Furstenau, Phys. Lett. B 260 (1991) 
447; P. Langacker and M.-x. Luo, Phys. Rev. D 44 (1991) 817; C. Giunti, C. W. Kim 
and U. W. Lee, Mod. Phys. Lett. A 6 (1991) 1745. 

[55] H. Murayama and A. Pierce, Phys. Rev. D 65 ( 2002 ) 055009 |hep-ph/0108104|. 

[56] S. A. R. Ellis and J. D. Wells, Phys. Rev. D 91 (2015) 7, 075016 [arXiv: 1502.01362 
[hep-ph]]. 

[57] J. Ellis, T. Falk, and K.A. Olive, Phys. Lett. B444 (1998) 367 |arXiv:hep-ph/9810360|; 

J. Ellis, T. Falk, K.A. Olive, and M. Srednicki, Astr. Part. Phys. 13 (2000) 181 
[Erratum-ibid. 15 (2001) 413] |arXiv:hep-ph/990548l|; R. Arnowitt, B. Dutta and 
Y. Santoso, Nucl. Phys. B 606 (2001) 59 |arXiv:hep-ph/0102181|; M. E. Gomez, 

G. Lazarides and C. Pallis, Phys. Rev. D D61 (2000) 123512 |arXiv:hep-ph/9907261|; 
Phys. Lett. B487 (2000) 313 |arXiv:hep-ph/0004028]; Nucl. Phys. B B638 (2002) 165 
|arXiv:hep-ph/0203131|; T. Nihei, L. Roszkowski and R. Ruiz de Austri, JHEP 0207 
(2002) 024 |arXiv:hep-ph/0206266]. 

[58] J. L. Feng, K. T. Matchev and T. Moroi, Phys. Rev. Lett. 84, 2322 (2000) |arXiv:hep- 
ph/9908309)]; Phys. Rev. D 61, 075005 (2000) |arXiv:hep-ph/9909334|; J. L. Feng, 

K. T. Matchev and F. Wilczek, Phys. Lett. B 482, 388 (2000) |arXiv:hep-ph/0004043|; 

H. Baer, T. Krupovnickas, S. Profumo and P. Ullio, JHEP 0510 (2005) 020 |hep- 
ph/0507282j; 

[59] C. Boehm, A. Djouadi and M. Drees, Phys. Rev. D 62, 035012 (2000) |arXiv:hep- 
ph/9911496||; J. R. Ellis, K. A. Olive and Y. Santoso, Astropart. Phys. 18, 395 (2003) 
|arXiv:hep-ph/0112113|; J. Edsjo, M. Schelke, P. Ullio and P. Gondolo, JCAP 0304, 001 
(2003) |arXiv:hep-ph/0301106|; J. L. Diaz-Cruz, J. R. Ellis, K. A. Olive and Y. Santoso, 
JHEP 0705, 003 (2007) |arXiv:hep-ph/0701229|; I. Gogoladze, S. Raza and Q. Shah, 
Phys. Lett. B 706, 345 ( 2012 ) |arXiv:1104.3566 [hep-ph]]; M. A. Ajaib, T. Li and 
Q. Shaft, Phys. Rev. D 85, 055021 ( 2012 ) [arXiv: 1111.4467 [hep-ph]]. 

[60] P. Skands et al. , JHEP 0407 (2004) 036 |arXiv:hep-ph/0311123|; B. Allanach et al, 
Comput. Phys. Commun. 180 (2009) 8 |arXiv:0801.0045 [hep-ph]]. 


47 



